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SUMMARY 


The  use  of  a  comprehensive  analytical  model  of  rotorcraft 
aerodynamics  and  dynamics  is  described.  This  analysis  is  designed 
to  calculate  rotor  performance,  loads,  and  noise;  the  helicopter 
vibration  and  gust  response;  the  flight  dynamics  and  handling 
qualities;  and  the  system  aeroelastic  stability.  The  analysis  is 
a  combination  of  structural,  inertial,  and  aerodynamic  models,  that 
.is  applicable  to  a  wide  range  of  problems  and  a  wide  class  of 
vehicles.  The  analysis  is  intended  for  use  in  the  design,  testing, 
and  evaluation  of  rotors  and  rotorcraft,  and  to  be  a  basis  for 
further  development  of  rotary  wing  theories.  This  report  describes 
the  use  of  the  computer  program  that  implements  the  analysis. 

/> 
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1.  PROGRAM  SUMMARY 


The  computer  program  calculates  the  loads  and  motion  of  helicopter 
rotors  and  airframe.  First  the  trim  solution  is  obtained;  then  the  flutter, 
flight  dynamics,  and/or  transient  behavior  can  be  calculated.  Either  a 
new  job  can  be  initiated,  or  further  calculations  can  be  performed  for  an 
old  job. 

For  a  new  job,  the  input  consists  of  block  data  or  an  input  file 
(the  program  can  create  the  input  file  from  the  block  data),  and  airfoil 
files.  Then  namelists  are  read  for  additional  data,  particularly  case- 
specific  inputs.  One  or  more  cases  can  be  run  for  a  new  job. 

For  an  old  job,  the  input  consists  of  a  restart  file  (written  during 
the  execution  of  a  previous  job) ,  and  name] ists .  Only  one  case  can  be  run 
for  an  old  job.  The  job  can  be  resumed  either  at  the  point  where  the  trim 
solution  vas  completed,  or  it  can  be  resumed  in  one  of  the  subsequent  tasks. 
For  a  trim  restart,  any  or  all  of  the  other  tasks  can  be  initiated.  For 
flutter,  flight  dynamics,  or  transient  restarts,  only  that  task  can  be  done. 

For  both  new  and  old  jobs,  a  scratch  file  is  usually  needed;  and  the 
job  may  write  data  on  the  restart  file.  In  the  flutter  and  flight  dynamics 
tasks,  eigenvalue  data  -,?y  be  written  on  a  file. 

For  both  new  and  old  jobs,  a  case  namelist  is  always  read  to 
define  the  job,  and  a  trim  namelist  is  read  to  define  the  flight  condition 
and  analysis  tasks.  Component  and  task  namelists  may  be  read  as  required. 

The  loads  and  motion  solution  is  obtained  by  an  iterative  process. 

The  inner-most  loop  consists  of  the  rotor  and  airframe  motion  calculation, 
for  prescribed  control  positions,  induced  velocity  distribution,  and  mean 
shaft  motion.  Convergence  of  the  motion  solution  is  determined  by  comparing 
the  calculated  harmonics  every  few  revolutions.  The  next  loop  consists  of 


the  uniform  or  nonuniform  rotor-induced  velocity  calculation,  followed 
by  the  motion  solution.  Convergence  is  determined  by  comparing  the  rotor 
th'-jst  or  circulation  used  to  calculate  the  induced  velocity  with  that 
resulting  after  the  motion  has  been  re-calculated.  Before  beginning  the 
circulation  and  motion  iterations,  the  blade  bending  and  torsion  modes  are 
calculated.  If  the  rotor  nonuniform  induced  velocity  is  used,  there  is 
an  additional  outer  loop,  consisting  of  calculation  of  the  rotor  wake 
influence  coefficients  followed  by  the  circulation  and  motion  iterations. 

To  calculate  the  influence  coefficients,  the  prescribed  or  free  wake  geometry 
must  be  evaluated.  Having  completed  the  motion  solution,  the  performance, 
loads,  vibration,  and  noise  can  be  evaluated  as  required. 

The  trim  analysis  proceeds  in  stages.  In  the  first  stage  the  trim 
solution  is  obtained  for  uniform  inflow;  in  the  second  and  third  stages  the 
trim  solution  is  obtained  for  nonuniform  inflow,  with  prescribed  or  free 
wake  geometry  respectively.  The  analysis  can  stop  at  any  of  these  stages. 

Within  each  stage,  the  aircraft  controls  and  orientation  are  incremented 
until  the  equilibrium  of  forces  required  for  the  specified  trim  state  is 
achieved. 

In  the  flutter  analysis,  the  matrices  are  constructed  that  describe 
the  linear  differential  equations  of  motion,  and  the  equations  are  analyzed. 
Optionally  the  equations  are  reduced  to  just  the  aircraft  rigid  body 
degrees  of  freedom  (by  a  quasistatic  reduction),  and  the  equations  are 
analyzed  as  for  the  flight  dynamics  task. 

In  the  flight  dynamics  analysis,  the  stability  derivatives  are 
calculated  and  the  matrices  are  constructed  that  describe  the  linear  differential 
equations  of  motion.  These  equations  are  analyzed  (optionally  including 
a  numerical  Integration  as  for  the  transient  analysis). 

In  the  transient  analysis,  the  rigid  body  equations  of  motion  are 
numerically  integrated,  for  a  prescribed  transient  gust  or  control  input. 


2.  SUBPROGRAM  FUNCTIONS 


The  following  pages  list  the  subprograms  that  constitute  the 
analysis,  and  state  the  primary  function  of  each  subprogram.  Only  the 
subprograms  for  rotor  #1  are  listed;  the  subprograms  for  rotor  #2  have 
identical  functions. 


Subprogram 

Name 

MAIN  Primary  job  and  analysis  control 

TIMER  Program  timer 

INPTN  Input  for  new  job 

INPTO  Input  for  old  job 

INPTA1  Read  airfoil  table  file 

INPTR1  Read  rotor  namelist 

INPIVl  Read  wake  namelist 

INFIB  Read  body  namelist 

INF  .Til  Read  loads  namelist 

INPTF  Read  flutter  namelist  for  new  job 

INPTS  Read  flight  dynamics  namelist  for  new  job 

INPTT  Read  transient  namelist  for  new  job 

INPTG  Read  flutter  namelist  for  old  job 

INPTU  Read  flight  dynamics  namelist  for  old  job 

INPTV  Read  transient  namelist  for  old  job 

FILEI  Read  or  write  input  file 

FILEJ  Read  or  write  trim  data  file 

FILER  Read  or  write  restart  file 

FILEF  Read  or  write  flutter  restart  file 

FILES  Read  or  write  flight  dynamics  restart  file 

FILET  Read  or  write  transient  restart  file 

FILEE  Write  eigenvalue  file 

INTT  Initialization 

INITA  Initialize  environment  parameters 

INITC  Initialize  case  parameters 

INITR1  Initialize  rotor  parameters 

INITB  Initialize  airframe  parameters 

INITE  Initialize  drive  train  parameters 

CHEKR1  Check  for  fatal  errors 
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Subprogram 

Name 


PRNTJ 

PRNTC 

PENT 

PRNTR1 

PRNW1 

PRNTB 

PRNTF 

PRNTS 

PRNTT 

PRNTC 

TRIM 

TRIMI 

TRIMP 

FLUT 

FLUTM 

FLUTB 

FLUTR1 

FLUTIi 

FLUTA1 

FLUTL 

STAB 

STABM 

STABD 

STABB 

STABL 

STABP 

TRAN 

TRANI 

TRANP 

TRANC 

CONTRL 

GUSTU 

GUSTO 

PBRF 

PERFR1 


Print  job  input  data 

Print  case  input  data 

Print  trim  input  data 

Print  rotor  input  data 

Print  wake  input  data 

Print  body  input  data 

Print  flutter  input  data 

Print  flight  dynamics  input  data 

Print  transient  input  data 

Print  transient  gust  and  control  input  data 

Trim 

Calculate  trim  solution  by  iteration 
Print  trim  solution 

Flutter 

Calculate  flutter  matrices 

Calculate  flutter  aircraft  matrices 

Calculate  flutter  rotor  matrices 

Calculate  flutter  Inertia  coefficients 

Calculate  flutter  aerodynamic  coefficients 

Analyze  flutter  constant  coefficient  linear  equations 

Flight  dynamics 

Calculate  flight  dynamics  stability  derivatives  and  matrices 

Print  stability  derivatives 

Calculate  flight  dynamics  equations 

Analyze  flight  dynamics  linear  equations 

Print  flight  dynamics  transient  solution 

Transient 

Calculate  transient  acceleration  for  numerical  integration 

Print  transient  solution 

Calculate  transient  gust  and  control 

Calculate  transient  control  time  history 

Calculate  uniform  gust  time  history 

Calculate  convected  gust  wave  shape 

Performance 

Calculate  and  print  rotor  performance 


Subprogram 

Name 


LOAD 

LOADR1 

L0ADH1 

L0ADS1 

LOADI1 

LOADF 

LOADH 

GEOMP1 

POLRPP 

HISTPP 

NOISR1 

BESSEL 


Loads,  vibration,  and  noise 

Calculate  and  print  rotor  loads 

Calculate  and  print  hub  and  control  loads 

Calculate  and  print  blade  section  loads 

Calculate  inertia  coefficients  for  section  loads 

Calculate  fatigue  damage 

Calculate  mean  and  half  peak -to- peak 

Printer- plot  of  wake  geometry 

Printer-plot  of  polar  plot 

Printer-plot  of  azimuthal  time  history 

Calculate  and  print  far  field  rotational  noise 

Calculate  J  Bessel  function 


RAMF 


Calculate  rotor/airframe  periodic  motion  and  forces 


M0DE1 

MO EEC 1 

MO  DEB  1 

MO  DEG 

MODEAl 

M0DET1 

MODEK1 

M0DED1 

INRTC1 

M0DEP1 


Blade  modes 

Initialize  blade  mode  parameters 

Calculate  blade  bending  modes 

Calculate  Galerkin  functions  for  bending  modes 

Calculate  articulated  blade  flap  and  lag  modes 

Calculate  blade  torsion  modes 

Calculate  kinematic  pitch-bending  coupling 

Calculate  blade  root  geometry 

Calculate  blade  inertia  coefficients 

Print  blade  modes 


BODYC 

ENGNC 

M0TNC1 

B0DYM1 

ENGNM1 


Initialize  airframe  parameters  at  trim 
Initialize  drive  train  parameters  at  trim 
Initialize  rotor  parameters  at  trim 
Calculate  airframe  transfer  function  matrix 
Calculate  drive  train  transfer  function  matrix 


WAKEU1 

WAKEN1 

INRTM1 

INRTI 

wnrmm* 
i  iv/  xiUiJ. 


Calculate  uniform  wake-induced  velocity 
Calculate  nonuniform  wake-induced  velocity 


Calculate  rotor  transfer  function  matrix 
Calculate  inverse  of  transfer  function  matrix 
Calculate  harmonics  of  hub  motion 


M0TNR1  Calculate  harmonics  of  rotor  motion 

M0TNB1  Calculate  blade  and  hub  motion 

AER0F1  Calculate  blade  aerodynamic  forces 

AER0S1  Calculate  blade  section  aerodynamic  coefficients 

AER0T1  Interpolate  airfoil  tables 


B0DYV1 

ENGNV1 

MOTNFi 

MOTNS 

BODYF 

BODYA 


Calculate  harmonics  of  airframe  motion 
Calculate  harmonics  of  drive  train  motion 
Calculate  rotor  generalized  forces 
Calculate  static  elastic  motion 
Calculate  airframe  generalized  forces 
Calculate  body  aerodynamic  forces 
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Subprogram 

Name 

WAKEC1  Calculate  influence  coefficients  for  nonuniform  inflow 

WAKE31  Calculate  blade  position 

VTXL  Calculate  vortex  line  segment  induced  velocity 

VTXS  Calculate  vortex  sheet  segment  induced  velocity 

GE0ME1  Evaluate  wake  geometry 

GE0MR1  Calculate  wake  geometry  distortion 

GE0MF1  Calculate  free  wake  geometry  distortion 

MINV  Calculate  inverse  of  matrix 

MINVC  Calculate  inverse  of  complex  matrix 

S1GENJ  Calculate  eigenvalues  and  eigenvectors  of  matrix 

DEREO  Eliminate  equations  and  variables  from  system  of  differential 

equations 

QSTRAN  Quasistatio  reduction  of  system  of  linear  differential  equations 

CSYSAN  Analyze  system  of  constant  coefficient  linear  differential 

equations 

DETRAN  Transform  equations  to  state  variable  form 

SINE  Calculate  frequency  response  from  matrices 

STATIC  Calculate  static  response  from  matrices 

ZERO  Calculate  zeros 

ZETRAN  Transform  matrix  for  zero  calculation 

BODE  Calculate  and  printer-plot  transfer  function  (Bode  plot) 

BODEPP  Printer-plot  transfer  function  magnitude  and  phase 

TRACKS  Calculate  and  printer- plot  time  history  of  time- invariant 
system  response 

TRCXPP  Printer-plot  time  history 

GUSTS  Calculate  and  print  rms  gust  response 

PSYSAN  Analyze  system  of  periodic  coefficient  linear  differential 
equations 

DEPRAN  Transform  equations  to  state  variable  form 
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3.  NAMELIST,  FILE,  AND  COMMON  BLOCK  LABELS 


The  list  below  gives  the  namelist  labels  used  by  the  program, 
and  the  type  of  input  data  read  in  each.  The  corresponding  common  block 
labels  are  given  in  the  right-hand  column. 


Namelist 

Common  Block 

Label 

Label 

NLCASE 

Job  data 

NLTRIM 

Trim  data 

TMDATA 

NLRTR 

Rotor  data 

R1DATA 

NLWAKE 

Wake  data 

G 1  DATA,  W1  DATA 

NLBODY 

Airframe  and  drive  train  data 

BDDATA  .BADATA  .ENDATA 

NLLOAD 

Loads  data 

LA DATA , LI DATA 

NLFLUT 

Flutter  data 

FLDATA 

NLSTAB 

Flight  dynamics  data 

STDATA,GCDATA 

NLTRAN 

Transient  data 

TNDATA  ,GCDATA 

The  list  below  gives  the  files  used  by  the  program.  The  left-hand 


column  gives 

the  input  parameter 

that  defines  the  file  unit  number. 

Unit 

Number 

NFDAT 

Input  data 

UFAFl 

Rotor  #1  airfoil 

data 

NFAF2 

Rotor  # 2  airfoil 

data 

Nl-RS 

Restart  data 

NFEIG 

Eigenvalue  data 

NFSCR 

Scratch  data 

The  list  below  gives  the  labels  of  the  common  blocks  used  by 
the  program,  and  states  the  type  of  data  contained  in  each.  Only  the 
common  blocks  for  rotor  #1  are  listed;  the  common  blocks  for  rotor  #2 
have  identical  functions. 
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Common  Block 
Label 


TMDATA 
R1DATA 
W1DATA 
G1DATA 
BDDATA 
BADATA 
ENDATA 
L1DATA 
LA  DATA 
GCDATA 
TNDATA 
STDATA 
FLDATA 
A1TABL 
UNITNO 
CASECM 

TRIMCM 

RTR1CM 

RH1CM 

BODYCM 

ERGNCM 

GUSTCM 

CONTCM 

CONVCM 

MD1CM 

INC1CM 

WKV1CM 

MNH1CM 

AES 1 CM 

MNR1CM 

MNSCM 

AEF1CM 

QR1CM 

r\T>  TV~»U 

i 

WG1CM 

WKC1CM 

AEMNCM 

LDMNCM 

FLMCM 

FLM1CM 

FLMACM 

FLINCM 

FLAECM 

STDCM 

STMCM 

TRANCM 


Input  trim  data 

Input  rotor  data 

Input  wake  data 

Input  free  wake  geometry  data 

Input  airframe  data 

Input  airframe  aerodynamics  data 

Input  drive  train  data 

Input  rotor  loads  data 

Input  airframe  loads  data 

Input  gust  and  control  data 

Input  transient  data 

Input  flight  dynamics  data 

Input  flutter  data 

Rotor  airfoil  tables 

Input/output  unit  numbers 

Job  description 

Calculated  trim  data 
Calculated  rotor  data 
Transfer  function  matrices 
Calculated  airframe  data 
Calculated  drive  train  data 
Gust  and  transient  control 
Aircraft  controls  and  motion 
Circulation  and  motion  convergence 
Blade  modes 

Rotor  inertial  coefficients 
Induced  velocity 
Hub  motion 

Blade  section  aerodynamics 
Rotor  motion  and  airframe  vibration 
Static  elastic  motion 
Rotor  forces 

Rotor  generalized  forces 

k  4  M  A  i-rii  nn-n'i  *1  A  J  .  |~1  r3 

nxi.ij.umo  goiioxUxx^ou  a \jj-  v/Gu 

Wake  geometry 

Wake  influence  coefficients 

Calculated  motion  data 
Calculated  loads  data 
Flutter  matrices 
Flutter  rotor  matrices 
Flutter  airframe  matrices 
Flutter  inertial  coefficients 
Flutter  aerodynamic  coefficients 
Flight  dynamics  stability  derivatives 
Flight  dynamics  matrices 
Calculated  transient  data 
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4.  PROGRAM  SKELETON 

The  following  pages  present  a  schematic  of  the  program,  showing 
the  basic  flow  of  control  and  the  major  loops,  options,  and  decisions. 

The  appearance  of  a  subprogram  name  (always  in  capital  letters)  means 
that  the  subprogram  is  called  at  that  point  in  the  analysis.  The  contents 
of  a  subprogram  are  shown  by  means  of  a  three-sided  box  appearing  below 
the  subprogram  name.  The  schematic  defines  the  input  and  output  structure 
of  the  program.  Timer  calls  and  trace-debug  prints  are  also  shown. 
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MAIN 

I  read  namelist  NLCASE 
if  new  job  and  BLKDAT  >  0 
DATE  (for  FILEID) 

TIME  (for  FILEID) 

FILEI  (input  file  write) 

PRNTJ 

for  JCASE  =  1  to  NGASES 

TIMER  (initialize) 

TIMER 

DATE  (for  I DENT) 

TIME  (for  IDENT) 
if  new  job 

INPTN 

INIT 

INITA 

INITC 

INITR1 

INITR2 

INITB 

INITE 

GHEKR1 

CHEKR2 

if  old  job 

INPTO 

PRNTG 

if  new  job  or  trim  restart 
TRIM 

FILEJ  (trim  data  scratch  file  write) 
if  ANTYPE(l)  /  0  or  flutter  restart 
FLUT 

FILEJ  (trim  data  scratch  file  read) 
if  ANTYPE(2)  /  0  or  flight  dynamics  restart 
STAB 

FILEJ  (trim  data  scratch  file  read) 
if  ANTYPE(3)  /  0  or  transient  restart 
TRAN 

TIMER 

TIMER  (print) 


INPTN 

FILEI  (input  file  read) 
read  namelist  NLTRIM 
if  OPREAD(l)  t  0 
INPTR1 


[read  namelist  NLRTR 

if  0PREAD(2)  t  0 
INPW1 


(read  namelist  NLWAKE 

if  0PREAD(3)  t  0 
INPTR2 


{read  namelist  NLRTR 

if  0PREAD(4)  /  0 
INPTW2 


jread  namelist  NLWAKE 

if  0PREAD(5)  t  0 
INPTB 


jread  namelist  NLBODY 

if  OPREAD(6)  *  0 
INPTL1 


[read  namelist  NLLOAD 

if  OHidAD(7)  /  0 
INPTL2 


jread  namelist  NLLOAD 

if  0PREAD(8)  t  0 
INPTF 


jread  namelist  NLFLUT 

if  OPREAD(9)  /  0 
INPTS 


jread  namelist  NLSTAB 

if  OPREAD(IO)  t  0 
INPTT 


read  namelist  NLTRAN 


if  first  case 

INPTA1 


read  airfoil  #1  file 


INPTA2 


BLKDAT.REFILE 


V 


...  - 


INPTO 

FILER  (restart  file  read) 

[filei 

FILEJ 

FILEF 

FILES 

FILET 


read  namelist  NLTRIM 
if  OPREAD(6)  }i  0 
INPTL1 


{read  namelist  NLLOAD 
if  OPREAD(7)  *  0 
INPTL2 


[read  namelist  NLLOAD 
if  OPREAE(S)  /  0  ~— 

INPTF 


read  namelist  NLFLUT 
INPTG 


[read  namelist  NLFLUT 
if  OPREAD(9)  *  0 
INPTS 

gT  namelist  NLSTAB 
INPTIJ 

read  namelist  NLSTAB  flutter  or 

if  OPREAD(lO)  /  0  . 

INPTT 


[read  namelist  NLTRAN 
INPTV 


read  namelist  NLTRAN 


flutter  restart 
flight  dynamics  restart 
transient  restart 


trim  restart 
flutter  restart 

trim  restart 
flight  dynamics  restart 

trim  restart 
transient  restart 
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TIMER 


NPRNTP  >  0 
NPRNTL  >  0 


LEVEL(l) 
LEVEL( 2) 


NPRNTP 

NPRNTL 


LEVEL( 1 ) 
LEVEL( 2) 


NPRNTP  >  0 
NPRNTL  >  0 

RSWRT  /  0 

IFILEI 
|  FILE J 


trim  restart  entry  point 
PR  NT 


PRNTC 

if  NFRNTI  t  0 
FRNTR1 
PRNTW1 
PRNTE2 
PRN1W2 
PRNTB 


M0EEP1 

M0DEP2 

TRIMP 

PERF 

LOAD 

TIMER 


if  trim  restart,  go  to  restart  entry  point 

uniform  inflow 
if  ITERU  0 

TRIMI 

if  NPRNTT  -  1 
PERF 
LOAD 

nonuniform  inflow,  prescribe  wake  geometry 
for  IT  =  1  to  ITERR 
WAKEC1 
WAKEC2 
TRIMI 

if  IT  =  multiple  NPRNTT 
PERF 
LOAD 

nonuniform  inflow,  free  wake  geometry 
for  IT  =  1  to  ITERF 
WAKEC1 
WAKEC2 
TRIMI 

if  IT  =  multiple  NPRNTT 
PERF 
LOAD 

FILE®  (restart  file  write) 


\V  W  y  V  W  \V 


TRIMI 


RAMF 

if  MTRIM  -9F  0  or  no  trim  iteration,  re 'urn 
if  DEBUG(4)  1,  print  trim  iteration 
for  COUNTT  =  1  to  MTRIM 

if  COUNTT- 1  =  multiple  MTi.IMD,  construct  D"1 
for  I  =  1  to  MT 

increment  controls 
RAMF 

MINV 

Increment  controls 
RAMF 

if  DEBUG(4)»£  1,  print  trim  iteration 
test  trim  convergence 


OPTRIM 

OPrRIM 

EPTRIM, OPTRIM 


PERF 

TIMER 

PERFR1 

PERFR2 

TIMER 


L0ADR1 

L0ADR2 

TIMER 


L 


LOADR1 


MOTNB1 

if  MALOAD  ^  0 
GEOME1 
HISTPP 
GEOMP1 
POLRPP 
HISTPP 

if  MHLOAD  0 
LOADH1 

LOADM 

LOADF 

HISTPP 


for  IR  =  1  to  MR LOAD 
LOADS1 


[LOADI1 

LOADM 

LOADF 

HISTPP 


for  IN  =  1  to  MNOISE 
NOISR1 


[BESSEL 


NPLOT( 1-4) 
MWKGMP 

nplot(5-6?) 

NPLOT(5-67) 


NPLOT( 68-71) 


NPLOT(72-75) 
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FLUT 

TIMER 

for  OPFLOW  0  (constant  coefficients) 

if  flutter  restart,  go  to  restart  entry  point 
FLUTM 

FILEF  (restart  file  write)  RSWRT  f  0 

flutter  restart  entry  point 

FRNTF 

M0DEP1 

M0DEP2 

FLUTL 


TIMER 

CSYSAN 

ANTYPE(l) 

*  0 

FILEE  (eigenvalue  file  write) 

BODE 

ANTYPE(2) 

t  0 

TRACKS 

ANTYPE(3) 

0 

GUSTS 

ANTYPE(4) 

t  0 

TIMER 

if  OPFDAN  f  0 

STABD 

STABE 


for  OPFLOW  >  0  (periodic  coefficients) 
for  NT  -  0  to  MPSIPC 
FLUTM 

if  NT  =  MPSIPC 
FRNTF 
M0DEP1 
M0DEP2 

PSYSAN 

if  NT  -  MPSIPC 

FILEE  (eigenvalue  file  write) 
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FLUTM 


MO  DEI 

M0IE2 

FLUTR1 

FLUTR2 

FLUTB 


BODYF 


FLUTR1 


NB  ■  NBLADB  if  OPFLOW  >0,  1  if  OPFLOW  -  0,  MPSICC  if  OPFLOW  < 
for  JPSI  -  1  to  NB 
FLUTI1 
FLUTA1 
1 


for  IR  “  1  to  MRA 
ABR0S1 


STAB 


TIMER 
PR  NTS 

if  flight  dynamics  restart,  go  to  restart  entry  point 
STABM _ 

for  ID  B  i  to  21 

increment  controls  or  motion 
for  IT  =  1  to  ITERS 


WaKECI  level( 1 ) 
UAKEC2  LEVEL(2) 
RAMF 

PERF  NPRNTP 
LOAD  NPRNTL 


FILES  (restart  file  write)  RSWRT  J  0 


flight  dynamics  restart  entry  point 

STABD 

STABE 

TIMER 


STABE 


for  IEQ  =  1  to  12 
DERED 
STABL 


TIMER 

CSYSAN 

FILEE  (eigenvalue  file  write) 

BODE 

TRACKS 

GUSTS 

numerical  integration 
HINV 
STABP 


PRNTC 


for  IT 


1  to  TMAX/TSTEP 
TRANC 


CONTRL 

GUSTO 

GUSTC 


if  IT  *  multiple  NPRNTT 
STABP 


TRCKPP 


TIMER 


EQTYPE(IEQ)  +  0 


ANTYPE(l)  /  0 

ANTYPE(2)  *  0 
ANTYPE( 3)  *  0 
ANTYPE(4)  ^  0 
ANTYPE(5)  j*  0 


OPTRAN  =  1 
OPTRAN  -  2 
OPTRAN  =  3 


19- 


V  V  aV  V 


if  transient  restart,  go  to  restart  entrv  point 

MINV 

TRANP 

for  IT  -  1  to  TMAX/TSTEP 
TRANC 

CONTRL 

GUSTU 

GUSTG 

TRANI 

for  IT  =  1  to  ITERT 
WAKECi 
WAKEC2 
RAMF 

if  IT  =  multiple  NPRNTT 
TRANP 
PERF 
LOAD 

if  IT  =  raul  » '.pie  NRSTRT 

FI7ET  (restart  file  write) 
transient  restart  entry  point 

TRCKPP 


OPTRAN  «  1 
OPTRAN  =  2 
OPTRAN  *  3 


LEVEL 

LEVEL' 


NPRNTP  >  0 
NPRHTL  >  0 

RSWRT  /  0 


131! 


RAMF 


TIMER 
BODYC 
MOTNC1 
MO  DEI 
BCDYM1 


INRTI 


M0TNC2 

M0DE2 

BODYM2 


INRTI 

for  COUNTC  =  1  to  ITERC  (circulation  iteration) 
WAKEU1 
WAKEN1 
WAKEU2 
WAKEN2 

for  COUNTM  =  1  to  ITERM  (motion  iteration) 
INRTM1 


INRTI 


INRTM2 


INRTI 


ENGNC 

ENGNM1 


INRTI 


ENGNM2 


INRTI 

for  JPSI  -  0  to  MREV  *  MPSI  by  MPSIR  (^  loop) 
M0TNH1 
M0TNR1 
M0TNK2 
M0TNR2 
B0DYV1 
ENGNV1 
MOTNF1 
B0DYV2 
ENGNV2 
M0TNF2 
MOTNS 

test  motion  convergence 
test  circulation  convergence 

BODYF 

Ibodya ” 


TIMER 


EPMOTN 

EPOIRC 
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WAKEB1 

GE0ME1 

for  I  »  1  to  MPSI  (*V  loop) 

WAXEB1 

WAKEB2 

for  M  =  1  to  NBLADE  (blade  loop) 

GE0ME1 

VTXL 

for  K  =  1  to  KFW  or  KDW  (A  loop) 


INFL0W(3)  = 


CSYSAN 


detran 

eigenj 

SINE 

STATIC 

ZERO 


ZETRAN 

EIGENJ 


BODE 


I®  TRAN 

eigenj 

ZERO 

[ZETRAN 

EIGENJ 


[BODEPP 


TRACKS 


[DEITAN 

eigenj 

MINVC 

TRCKPP 


GUSTS 


DETRAN 

EIGENJ 

MINVC 


PSYSAN 


depran 

eigenj 
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5 .  JOB  STRUCTURE 


In  this  section  the  structure  of  a  job  to  run  the  program  is 
defined.  The  basic  structure  consists  of  the  following  steps: 

1)  File  definition  as  required  for  job 

2)  Block  data  load  for  airframe  and  each  rotor 

3)  Main  program  call 

4)  Namelist  &NLCASE 

5)  Namelist  &NLTRIM  (for  each  case) 

6)  Component  and  task  namelists  as  required 
File  definition  parameters: 

a)  RET  =  T  Erase  file  at  logoff 

b)  DISP  =  NEW  New  file  to  be  created 

c)  DISP  =  OLD  Existing  file 

Sample  jobs  are  presented  below. 


New  job,  2  cases s  trim  analysis?  block  data  input,  basic  namelist  input, 
same  airfoil  table  for  both  rotors 

DDEF  FT50F001 , , SCRATCH , DISP=NEW , RET=T 
DDEF  FT41F001 .AIRFOIL, DISP=OLD 
LOAD  HELAj  LOAD  HELR1;  LOAD  HELR2 
CALL  MAINPROG 

&NLCASE  J0B=0,NCASES”2,RSWRT=0,BLKDAT=-1 , 

NFAF1-41 , NFAF2=41 , NFSCR=50 , NFRS=-1 , NFEIG=- 1 , 

&END 

&NLTRIM  VKT5=x , ,COLL=x. ,LATCYC=x. ,LNGCYC=x. ,PEDAL=x. ,APITCH=x. ,AR0LL=x. , 
ANTYFE=3*0,OPREAD=iO*0, 

&END 

&NLTRIM  data  for  second  case ,&END 
#END 
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s* 


New  job,  1  case;  trim,  flutter,  flight  dynamics,  and  transient  analysis; 
block  data  input,  all  namelist  inputs,  different  airfoil  table  for  each 
rotor;  write  eigenvalue  file 

DDEF  FT50F001, .SCRATCH, DISP= NEW, RET=T 
DDEF  FT41F001 , ,  AIRF0IL1 ,  DI3POLD 
DDEF  FT42F001 , , AIRF0IL2, DISP=OLD 
DEEF  FT45F001 , , EIGEN, DISP=NEW 
LOAD  KELA;  LOAD  HELR1 ;  LOAD  HELR2 
CALL  MA INFROG 

&NLCASE  JOB=0,NCASES=1,RSWRT=0,BLKDAT~-1 , 

NFAF1=41 , NFAF2=42 , NFSCR=50 , NFRS=-1 ,NFEIG=45 , 

&END 

&NLTRIM  VKTS=x. , 

COLL^x. ,LATCYC=x. ,LNGCYC=x. ,PEDAL=x. ,APITCH=x. ,AROLL=x. , 
ANTYPE=3*1 ,0FREAD*10*1 , 

&END 

&NLRTR  data.&END 
&NLWAKE  data ,&END 
&NLRTR  data.&END 
&NLWAKE  data.&END 
&NLBODY  data.&END 
&NLLOAD  data.&END 
4NLL0AD  data.&END 
&NLFLUT  data.&END 
&NLSTAB  data.&END 
&NLTRAN  data.&END 
$END 


New  job,  1  case;  trim  analysis;  block  data  input  and  write  input  file 

DDEF  FT50F001,, SCRATCH, DISP=*NEW,RET=T 
DDEF  FT41F001 , .AIRFOIL , DISP=0LD 
DDEF  FT40F001 . . INPUT , DISP=NEW 
LOAD  HELA;  LOAD  HELR1 ;  LOAD  KELR2 
CALL  HAINFROG 

&NLCASE  J0B=0, NCASES=1 ,RSWRT=0,BLKDAT=1 , 

NFAFl=4l  ,NFAF2°4l  ,NFSCR=50,NFRS=-1  ,NFEIG=-1  .NFDAT^O, 

&END 

&NLTRIM  data,4END 
$END 
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New  job,  1  case;  trim  analysis;  read  input  file 

DDEF  FT50F001, .SCRATCH, DISP=NEW,RET=T 
DDEF  FT41F001 , , AIRFOIL , DiSP=OLD 
DDEF  FT40F001 , , INPUT, DISP=0LD 
CALL  MAINPROG 

&NLCASE  J0B==0,NCASES=1,RSWRT=0,BLKDAT=0,RDFILE=1, 
NFAFl=4l ,NFAF2=4l ,NFSCR=50, NFRS=-1 ,NFEIG=-1 ,NFDAT=40, 
&END 

&NLTRIM  data,&END 
^END 


New  job,  2  cases;  trim  and  flutter  analysis;  write  restart  file 

DDEF  FT50F001 , , SCRATCH , DISP=NEW , RET=T 
DDEF  FT41F001, .AIRFOIL, DISP=0LD 
DDEF  FT44F001 ,  .RESTART1 ,  DISP=NEW 
DDEF  FT44F002,,RESTART2,DISP=NEW 
LOAD  HELA;  LOAD  HELR1 ;  LOAD  HELR2 
CALL  MAINPROG 

&NLCASE  J0B= 0 , NCA SES= 2 , RS WR 1 ,BLKDAT=-1 , 

NFAFl=4l .NFAF2-41 ,NFSCR=50,NFEIG=-1 , NFRS=44, 

*END 

&NLTRIM  data  for  first  case, 

ANTYPE=1 , 0 , 0 ,OPREAD( 8)=1 , 

&END 

&NLFLUT  data,&END 
&NLTRIM  data  for  second  case,&END 
&NLFLUT  dataf&END 
%EM) 


Old  job;  trim  restart  with  flutter  analysis 

DDEF  FT50F001 , .SCRATCH, DIS?=NEW,RET=T 
DDEF  FT44F001,, RESTART, DTSP=0LD 
CALL  MAINPROG 

&NLCASE  J0B=1 ,RSWRT=1 ,START=1 , 
NFSCR*50,NFEIG=-1 ,NFRS=44, 

&END 

&NLTRIM  ANTYPE=l,0,0,OPREAD(8)=i, 
&END 

&NLFLUT  data,&END 
$END 
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Old  job;  flutter 


restart 


!}JJE  n?,0F001  ’  »scratch,  DISP=*NEW  .RET^T 

r\m  IrJ^F°01  ’  *RESTART,  DISP=OLD 
CALL  MA INPROG 

&NLCASE  JOB=l  ,RSWRT=0,  STARTS 
NFSCR=50 , NFEIG=- 1 , NFRS=44 
&END 

&NLTRIM  OPREAD(8)*=1 
&END  ’ 


&NLFLUT  data.&END 
%EKD 
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6.  INPUT  DESCRIPTION 

In  this  section  the  input  variables  for  the  program  are  defined. 
The  variables  are  categorized  according  to  the  namelist  that  reads  them. 
The  program  namelist  labels  are  listed  in  the  table  below. 


Namelist 

Label 

NLCASE 

Job  data 

NLTRIM 

Trim  data 

NLRTR 

Rotor  data 

NLWAKE 

Wake  data 

NLBODY 

Airframe  and  drive  train  data 

NLLOAD 

Loads  data 

NLFLUT 

Flutter  data 

NLSTAB 

Flight  dynamics  data 

NL1RAN 

Transient  data 

The  corresponding  common  block  labels,  for  the  block  data  form  of  input, 
may  be  obtained  from  Section  3*  In  the  description  of  the  input  parameters 
for  the  rotor,  the  variables  NBM  and  NTM  are  used: 

a)  NBM  is  the  index  of  the  highest-frequency  blade  bending 
mode  used  in  the  analysis; 

b)  NTM  is  the  index  of  the  highest-frequency  blade  torsion 
mode  used  in  the  analysis. 
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Namelist  NLCASE 


JOB 

RSWRT 


NCASES 

BLKDAT 


RDFILE 


START 


NFDAT 

NFAF1 

NFAF2 

NFRS 

NFEIG 

NFSCR 

NUIN 

NUOUT 

NUEB 

NUPP 

NULIN 


integer  parameter  defining  job:  EQ  0  for  new  job  (default); 
NE  0  for  old  job  or  restart  (one  case  only) 

integer  parameter  controlling  restart  file  write:  0  to 
suppress  write  (default) 


New  job  only 

number  of  cases  (default  =  l) 

integer  parameter  defining  input  source: 

EQ  0  read  input  file  (default) 

GT  0  use  loaded  block  data  and  write  input  file 

LT  0  use  loaded  block  data 

integer  parameter  controlling  input  file  read: 

EQ  0  read  file  for  first  case  only 

NE  0  read  file  for  every  case  (default) 

Old  job  only 

integer  parameter  defining  task: 

1  for  trim  restart  (default) 

2  for  flutter  restart 

3  for  flight  dynamics  restart 

4  for  transient  restart 

trim  restart  can  be  followed  by  any  or  all  of  the  other  tasks 
(as  defined  by  ANTYPE);  for  flutter,  flight  dynamics,  or 
transient  restart,  only  that  task  can  be  done 


Input/output  unit  numbers 

input  data  file  (new  job  only);  default  =  40 

rotor  #1  airfoil  file  (new  job  only);  default  =  41 

rotor  #2  airfoil  file  (new  job  only;  only  if  have  two  rotors); 
default  =  42 

restart  file  (no  file  write  if  LE  0);  default  =  44 

eigenvalue  file  (no  file  write  if  LE  0);  default  =  45 

scratch  file;  default  =  50 

namelist  input;  default  =  5 

printer  (and  debug  level  l);  default  =  6 

debug  output  (levels  2  and  3)i  default  =  6 

printer- plots;  default  =  6 

linear  system  analysis;  default  =  6 
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Namelist  NLTRIM 


OFREAD(lO)  integer  vector  defining  namelist  read  structure;  EQ  0 

to  suppress  read : 

components  (new  job  only) 

(1)  NLRTR,  rotor  #1 

(2)  IfLWAKE ,  rotor  4 1 

(3)  NLRTR,  rotor  til 

(4)  NLWAKE,  rotor  #2 

(5)  NLBODY 

tasks 

i.6)  NLLCAD,  rotor  ti  1 

(7)  NLLOAD,  rotor  #2 

(8)  NLFLUT 

(9)  NLSTAB 

(10)  NLTRAN 


NFRNTI 

ANTYPE(3) 


TITLE (20) 
CODE 

OPUNIT 

NROTOR 


Integer  parameter  controlling  input  data  print:  EQ  0 
for  short  form  print  only 

integer  vector  defining  tasks  for  new  job  or  trim 
restart;  EQ  0  to  suppress: 

( 1 )  flutter 

(2)  flight  dynamics 

(3)  transient 

title  for  job  and  case  (80  characters) 

alphanumeric  code  for  job  and  case  identification; 

4  characters 

integer  parameter  designating  unit  system;  1  for 
English  units  (ft-slug-sec) ;  2  for  metric  units  (m-kg-sec) 

number  of  rotors 
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NLTRIM 


DEBUG(25)  integer  vector  controlling  debug  printt 

0  no  debug  print 

1  trace  print 

2  low  level  print. 

3  high  level  print 

(1)  time  (sec)  at  which  debug  print  enabled 

(2)  input,  2-3  ( INPTx) 

(3)  initialization,  2  (INITG,  INITR,  INITB,  INITE) 

(4)  trim  iteration,  1-2  (TRIHl) 

(5)  loads,  2  (LOADI) 

(6)  flutter  matrices,  2-3  ( FLUTM) 

(?)  flutter  coefficients,  2-3  (FLUTI,  FLUTA) 

(8)  flight  dynamics,  2-3  ( STABM ,  STABE) 

(9)  transient,  2  (TRANl) 

(10)  rotor/airframe  motion  and  forces,  2-3  (RAMF) 

(11)  blade  modes,  2  (MODE,  MODEx) 

(12)  inertia  coefficients,  2  (INRTC) 

(13)  airframe  constants  and  matrices,  2  (BODYC,  ENGNC, 

MOTNC ,  BODYM,  ENGNM) 

(lM  induced  velocity,  2  (WAKEU,  WAKEN) 

(15)  rotor  matrices,  2-3  (INRTM) 

(16)  hub/airframe  motion  and  generalized  forces,  2 
(MOTNH,  BODYV,  ENGNV,  MOTNF,  MOTNS) 

(17)  rotor  motion,  2-3  (MOTNR) 

(18)  rotor  aerodynamics,  2-3  (AEROF) 

(19)  blade  section  aerodynamics,  3  (AEROS) 

(20)  body  forces  and  aerodynamics,  2  (BODYF) 

(21)  wake  influence  coefficients,  2  (WAKEC) 

(22)  vortex  line  and  sheet,  3  (VTXL,  VTXS) 

(23)  prescribed  wake  geometry,  2-3  (GEOMR) 

(24)  free  wake  geometry,  1-3  (GEOMF) 

(25)  timer,  1  (TIMER) 


VKTS 

VEL 


TIP 

RFM 


CFPENS 


ALTMSL 

TEMP 

DENSE 

OPGRND 

HAGL 

OPENGN 


NLTRIM 

aircraft  speed  V  (knots) 
velocity  ratio  ”//lR 

input  either  VEL  or  VKTS  by  namelist;  if  neither 
parameter  is  defined,  V  =  0  is  used 

rotor  tip  speed  i~2.  3  (ft/sec  or  m/sec) 

rotor  rotational  speed  (rpm) 

input  either  VTIP  or  RPM  by  namelist;  if  neither 
parameter  is  defined,  the  normal  tip  speed  VTIPN 
is  used:  rotor  #2  speed  is  calculated  from  the 
i-ear  ratio  TRATIC 

integer  parameter  defining  specification  of  aerodynamic 
environment:  if  1,  given  altitude  and  standard  day; 
if  2,  given  altitude  and  temperature;  if  3  given  density 
and  temperature 

altitude  above  mean  sea  level  (ft  or  m) ,  for  CPDENS  =  1  or  2 

air  temperature  (°F  or  °C),  for  OPDENS  =  ?.  or  3 

air  density  ^slug/ft^  or  kg/rn^),  for  OPDENS  =  3 

integer  parameter  controlling  ground  effect  analysis: 

EQ  0  for  out  of  ground  effect,  NE  0  for  in  ground  effect 

altitude  helicopter  center  of  gravity  above  ground  for 
ground  effect  analysis  (ft  or  m) 

integer  parameter  specifying  engine  state ;  1  for  autorotation 

(engine  inertia,  engine  damping,  and  throttle  control  torque 
zero;  no  engine  speed  degree  of  freedom);  2  for  engine  out 
(engine  damping  and  throttle  control  torque  zero);  0  for 
normal  operation 

wing  flap  angle  &F  (deg) 

for  free  flight,  trim  turn  rate  Vp  (deg/sec),  positive  to  right 


AFLAP 

RTURN 


NLTRIM 


COLL 

LATCYC 

LNGCYC 

PEDAL 

APITCH 

AROLL 


ACLIMB 

AYAW 


MPSI 

MPSIR 

MREV 

ITERH 

EPMOTN 

ITERC 

EPCIRC 


initial  values  of  controls  (trimmed  as  appropriate) 

collective  stick  displacement  £  or  <5^0  (deg), 
positive  up  0  govr 

lateral  cyclic  stick  displacement  %  (deg),  positive  left 

longitudinal  cyclic  stick  displacement  C  (deg),  positive 
aft  s 


pedal  displacement  (deg),  positive  to  righr 

for  free  flight,  aircraft  pitch  angle  (deg),  positive 
nose  up;  for  wind  tunnel,  rotor  shaft  angle  of  attack  0_,, 
(deg),  positive  nose  up  1 


for  free  flight,  aircraft  roll  angle  (deg),  positive 

to  right 


(  and  define  orientation  of  body  axes  relative  to 
earth  axes) 


for  free  flight,  aircraft  climb  angle  ©Fp  (deg),  positive  up 

for  free  flight,  aircraft  yaw  angle  H'pp  (deg),  positive  to 
right;  for  wind  tunnel,  test  module  yaw  angle  y  (deg), 
positive  to  right 

( 0p  and  S*  define  orientation  of  velocity  axes  relative 
to  eArth  axes;  Vclimb  =  V  sin  ©Fp  and  Vside  »  V sin*VFp  cos0pp) 


number  of  azimuth  steps  per  revolution  in  motion  and  loads 
analysis,  maximum  36;  for  nonuniform  inflow  must  be  multiple 
of  number  of  blades;  for  free  wake  geometry,  maximum  2*1- 

in  harmonic  motion  solution,  number  of  azimuth  steps  between 
update  of  airframe  vibration  and  rotor  matrices 

in  harmonic  motion  solution,  number  of  revolutions  between 
tests  for  motion  convergence 

maximum  number  of  motion  iterations 

tolerance  for  motion  convergence  (deg) 

maximum  number  of  circulation  iterations 

tolerance  for  circulation  convergence  (  /^C_/v) 


D0F(54) 


NLTRIM 


D0FT(8) 


MHARM(2) 


MHARMP(2) 


LEVEL ( 2) 


integer  vector  defining  degrees  of  freedom  used  in 
vibratory  motion  solution,  0  if  not  used;  order: 


rotor  #1 
rotor  ft 2 

airframe 


<1}  ...  q 


10 

l10 


?4 


(bending,  max  10)  (torsion,  max  5)  (gimbal/teeter) 


£ 


F  %  VF  *F  yF 
(rigid  body) 


drive  train 


(rotor/engine  speed) 


SF  • • •  4si6 

(flexible  body,  max  10) 

Afe  A©  A>© 

x.  govri  govr2 

(governor) 


integer  vector  defining  blade  bending  degrees  of  freedom 
usad  for  mean  deflection  (subset  of  DOF),  0  if  not  used; 
order: 

rotor  tt\  qt  q2  q^  q^ 

rotor  tt2  q1  qg  q^  q4 

(bending,  max  4) 


number  of  harmonics  in  rotor  motion  analysis;  maximum  20; 
EQ  0  for  mean  only 

(1)  rotor  #1 

(2)  rotor  itZ 

number  of  harmonics  in  airframe  vibration  analysis 
(harmonics  of  N/rev);  maximum  10;  EQ  0  for  static  elastic 
only;  suggest  LE  MHARH/NB4ADE,  and  the  same  value  for 
both  rotors  if  coupled  hub  vibration  used  (see  OFrtVIB) 

(1)  rotor  #1 

(2)  rotor  ttZ 

integer  parameter  specifying  rotor  wake  analysis  level: 

0  for  uniform  inflow,  1  for  nonuniform  inflow  with 
prescribed  wake  geometry,  2  for  nonuniform  inflow  with 
free  wake  geometry  (must  be  consistent  with  INFLOW) 

(1)  rotor  #1 

(2)  rotor  ttZ 
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ITERU 

ITERR 

ITERF 

NPRNTT 

NPRNTP 

NPRNTL 

MTRIM 

MTRIMD 

DELTA 

FACTOR 

EPTRIH 

OPGOVT 


CXTRIM 

XTRIM 

CTTR1M 

CPTRIM 

CYTRIM 

BCTRIM 

BSTRIM 


number  of  wake  and  trim  iterations 

at  uniform  inflow  levels  EQ  0  to  skip 

at  nonuniform  inflow/prescribed  wake  geometry  level: 

EQ  0  to  skip 

at  nonuniform  inflow/free  wake  geometry  level 

integer  parameter  ns  trim/performance/load  print 
every  n-th  iterations  LE  0  to  suppress 

integer  parameter  controlling  performance  print;  LE  0  to 
suppress 

integer  parameter  controlling  loads  prints  LE  0  to  suppress 

maximum  number  of  iterations  on  controls  to  achieve  trim 

number  of  trim  iterations  between  update  of  trim  derivative 
matrix 

contro!  step  in  trim  derivative  calculation  (stick  displacement, 

factor  reducing  control  increment  in  order  to  improve  trim 
convergence  (typically  0.5) 

tolerance  on  trim  convergence 

integer  parameter  specifying  governor  trim 
0  trim  collective  stick  & 

1  trim  rotor  #1  governor  0 

2  trim  rotor  #2  governor 

3  trim  both  rotor  governors 

targets  for  wind  tunnel  trim  cases 

v- 

X/q  (ft^  or  m^) 

cT/r 

Cp/<r 
Cy/o- 
£e  (deg) 

(deg) 


or  CL/v~ 


NLTR 


OPTRIM 


integer  parameter  specifying  trim  option 
free  flight  cases 


OFIRIM  =  0 
1 


a 


6 


no  trim 

trim  forces  and  moments  with 

K  &c  Ss  Sp  9  FT  ^FT 
trim  forces  and  moments  with 
So  Sc  &s  Sp  ®{TT  H>pp 
trim  forces,  moments,  and  power  with 

So  Sc  Ss  £> p  4>pT  &pp 

trim  forces,  moments,  and  power  with 
So  S3  &p  Qp*p  vp  ^pp 
trim  symmetric  forces  and  moments  with 
So  ^3  ^FT 

trim  symmetric  forces,  moments,  and  power  with 


So  Sg  Qp'p  fipp 


wind  tunnel 

cases 

cptr:m  =  10 

no  trim 

U 

trim 

cA 

with 

*o 

12 

trim 

C/r 

with 

eT 

13 

trim 

Cp/<r 

with 

*0 

14 

trim 

Pc  < 

2s 

with 

Sc 

Ss 

15 

trim 

<A 

Pc  ( 

2 

s 

with 

So 

Sc 

Ss 

16 

trim 

Vv 

ca 

cY/V 

with 

So 

Sc 

Ss 

17 

trim 

°iA- 

V0" 

cA 

with 

So 

Sc 

18 

trim 

ciA- 

G.A 

?c  P£ 

.  with 
> 

So 

Sc 

Ss  °T 

19 

trim 

ciA 

X/q  C 

Y/o" 

with 

So 

Sc 

Ss 

20 

trim 

CA 

X/q  C 

v/*“ 

JU 

with 

So 

Sc 

&r 

21 

trim 

X/q  | 

Fe  Ps 

with 

So 

Sc 

Ss 

22 

trim 

Pc 

with 

Ss 

23 

trim 

cA 

Pc 

with 

So 

Ss 

24 

trim 

cA 

ca 

with 

So 

Ss 

25 

trim 

Cl/r 

with 

So 

26 

trim 

cL/°" 

Pc 

with 

So 

Ss 

27 

trim 

Cj/c- 

X/q 

with 

So 

Ss 

28 

trim 

CL/<r 

X/q 

with 

Orp 

29 

trim 

CA 

X/q  j 

?c 

with 

So 

Ss 

©T 

-37- 


NLT°T^ 


WEIGHT  see  namelist  NLBCDY 
I XX 
IYY 
IZZ 
IX  Y 
IXZ 
IYZ 
ATILT 
FSCG 
BLCG 
WLCG 


-38- 


Namelist  NLRTR 


TITLE(20)  title  for  rotor  and  wake  data  (80  characters) 

TYPE  rotor  identification  (4  characters);  suggest  MAIN,  FRNT, 

or  RGHT  for  rotor  #1;  and  TAIL,  REAR,  or  LEFT  for  rotor  #2 

VTIPN  normal  tip  speed  nRQ  (ft/sec  or  m/sec) 

RADIUS  blade  radius  R  (ft  or  m) 

SIGMA  solidity  ratio  =  Nc  /'trR  (based  on  mean  chord) 

m  ji 

GAMMA  blade  Lock  number  ^  ^  ac  R  /i^  (based  on  standard 

density,  a  =  5*7,  ar.S  mean0 chord) 

( &  and  <r*  are  only  used  to  calculate  the  normalization 
parameters  c^  and  1^) 

NBLADE  number  of  blades 


TDAMPO 

TDAMPC 

TDAMPR 

NUGG 

NUGS 

GDAMPC 

GDAMPS 

LDAMPC 


T  rvAMTSM 

**  t*t» 

LDAMPR 

GSB(N3M) 

GST(NlM) 

ROTATE 


control  system  damping  (ft- lb/rad/sec  or  m-N/rad/sec) 
collective 


cyclic 


rotating 

longitudinal  gimbal  natural  frequency  ^G„  or  teeter  natural 
frequency  (per  rev  at  normal  tip  speed  VTIPN) 

lateral  gimbal  natural  frequency  ^GS  (per  rev  at  normal  tip 
speed  VTIPN) 

longitudinal  gimbal  damping  C„G  or  teeter  damping  C^, 

(ft- lb/rad/sec  or  m-N/rad/sec; 

lateral  gimbal  damping  CGS  (ft-lb/rad/sec  or  m-N/rad/sec) 

linear  lag  damper  coefficient  (ft-lb/rad/sec  or 
m-N/rad/sec);  estimated  damping  if  a  nonlinear  damper  is 
used  (LDAMPM  GT  0.);  the  lag  mode  has  structural  damping 
also  (GSB) 


maximum  moment  of  nonlinear  lag  damper-  M_ r  (ft- lb  or  m-N); 
linear  lag  damper  used  if  LDAMPM  EQ  0.  ^ 


lag  velocity  $T-  where  maximum  moment  of  lag  damper  occurs 
(rad/sec);  hydraulic  damping  below  and-  friction  damping 
above 


bending  mode  structural  damping  gg 
torsion  mode  structural  damping  g0 

o 

integer  parameter  specifying  rotor  rotation  direction: 
1  for  counter-clockwise,  -1  for  clockwise  (viewed  from 
above) 


NLRTR 


0PHVIB(3) 

BTIP 

OPTIP 

LINTW 

TWISTL 

OPUSLD 

OPCOMP 

INFLOW (6) 


RRCOT 

RGMAX 


integer  parameter  controlling  hub  vibration  contributions; 
gravity  and  static  velocity  terms  always  retained;  0  to 
suppress : 

(l)  vibration  due  to  this  rotor 
(?)  vibration  due  to  other  rotor  (must 
supress  if  Qg/fV  t  l) 

(3)  static  elastic  motion 

tip  loss  parameter  B 

integer  parameter  specifying  tip  loss  type:  1  for  tip  loss 
factor,  2  for  Prandtl  function 

integer  parameter  specifying  twist  type;  EQ  0  for  nonlinear 
twist,  HE  0  for  linear  twist 

linear  twist  rate  0,  (deg);  used  to  calculate  TWISTA  and 
WISTI  if  LINTW  NE  0 

integer  parameter  controlling  use  of  unsteady  lift,  moment, 
and  circulation  terms:  if  0,  suppress;  if  i,  include;  if  2, 
zero  for  stall  (15°  <  |  cx  l  <  165°) 

integer  parameter  controlling  aerodynamic  model,  EQ  0  for 
incompressible  loads 

Inflow  model 

integer  vector  defining  induced  velocity  calculation 
(must  be  consistent  with  LEVEL) 

(1)  at  this  rotor:  0  for  uniform,  1  for  nonuniform 

(2)  at  other  rotor:  0  for  zero,  1  for  empirical, 

2  for  average  at  hub,  3  for  nonuniform  (only 

if  xXj/rij  =  1) 

(3)  at  wing-body:  0  for  zero,  1  for  empirical, 

2  for  nonuniform 

(4)  at  horizontal  tail:  0  for  zero,  1  for 
empirical,  2  for  nonuniform 

(5)  at  vertical  tail:  0  for  zero,  1  for  empirical, 
2  for  nonuniform 

(6)  at  point  off  rotor  disk:  0  for  zero,  1  for 
nonuniform 

root  vortex  position  for  wake  model,  rroo^/^ 

rGmax^  (^-nc*ucec*  velocity  calculated  using  maximum  bound 
circulation  magnitude  outboard  of  rGmax) 
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HR  A 

RAE(MRA+  l) 


CHORD(MRA) 

XA(MRA) 

THETZL(MRA) 


twista(mra) 

xac(mra) 

mcorrl(mra) 

mcorrd(mra) 

mgorrh(mra) 


Blade  section  aerodynamic  characteristics 

number  of  aerodynamic  segments;  maximum  30 

radial  stations  r/R  at  edges  of  aerodynamic  segments; 
sequential,  from  root  to  tip 

Following  quantities  are  specified  at  midpoint  of 
aerodynamic  segment 

blade  chord,  c/R 

offset  of  aerodynamic  center  aft  of  elastic  axis,  xa/R; 
x.  is  the  point  about  which  the  moment  data  in  the 
airfoil  tables  is  given 

incremental  pitch  of  zero  lift  line,  O  (deg);  can  be 
included  in  TWISTA;  is  the  pitch  of  the  axis 

corresponding  to  zero  angle  of  attack  in  the  airfoil 
tables,  relative  to  the  twist  angle  (TVISTA) 

blade  twist  relative  .?5R*  (deg) 

offset  of  aerodynamic  center  (for  unsteady  aerodynamics) 
aft  of  elastic  axis,  x.  J'K 

AU 

Mach  number  correction  factor  f=M  M  for  lift 

M  ell 

Mach  number  correction  factor  f^  =  for  drag 

Mach  number  correction  factor  f „  =  M 


eff/M  for  moment 


MRI 

RI(MRI) 

MASS(MRI) 

EIXX(MRI) 

EIZZ(MRl) 

Xl(MRl) 

xc(mri) 

KP2(MRI) 

itheta(mri) 


Blade  section  inertial  and  structural  characteristics 

number  of  radial  stations  where  characteristics  defined; 
maximum  51 

radial  stations  r/R;  sequential,  from  root  to  tip, 

Rl(l)  =  0.  and  Rl(MRl)  =  1. 

section  mass,  m  (slug/ft  or  kg/m) 

chordwise  bending  stiffness  (lb-ft^  or  N-m^) 

flapwlse  bending  stiffness  (lb-fl^  or  N-m^) 

offset  of  center  of  gravity  aft  of  elastic  axis,  x^/R 

offset  of  tension  center  aft  of  elastic  axis,  Xp/R 
(at  the  tip,  XG  should  be  set  nearly  equal  XI) 

p  O 

polar  radius  of  gyration  about  elastic  axis,  kp/ R4" 

section  moment  of  inertia  about  elastic  axis,  1^ 
(slug-ft  or  kg-m) 

torsional  stiffness,  GJ  (lb-ft^  or  N-m^) 
blade  twist  relative  .75Ri  (deg) 
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GJ(MRT) 

TOISTI(MRI) 


NLRTR 


OPSTLL 


OPYAW 


TAU(3) 

A DELAY 

AMAXNS 

PSIDS(3) 

ALFDS(3) 

ALFRE(3) 

CLDSP 

CDDSP 

CMDSP 


Stall  model 

integer  parameter  defining  stall  model 
0  no  stall 

1  static  stall 

2  McCroskey  stall  delay 

3  McCroskey  stall  delay  with  dynamic  stall 
vortex  loads 

4  Boeing  stall  delay 

5  Boeing  stall  delay  with  dynamic  stall 

vortex  loads 

(the  stall  delay  can  be  suppressed  by  setting  TAU=  0.) 

integer  parameter  defining  yawed  flow  corrections 

0  both  yawed  flow  and  radial  drag  included 

1  no  yawed  flow  (cos A.  =  1.) 

2  no  radial  drag  (F  =0.) 

3  neither  yawed  flow  nor  radial  drag  included 

stall  delay  time  constants  for  lift,  drag,  and  moment: 

•C-,  n:D,  (calculated  if  LT  0.) 

maximum  angle  of  attack  increment  due  to  stall  delay, 

<x.  ,  (deg) 

max  delay  &/ 

angle  of  attack  in  linear  range  for  no  stall  model,  c><  (deg) 

^  max 

dynamic  stall  vortex  load  rise  and  fall  time  (azimuth  increment) 
for  lift,  drag,  and  moment:  .  (deg) 

CIS 

dynamic  stall  angle  of  attack  for  lift,  drag,  and  moment: 

^ds  (deg) 

stall  recovery  angle  of  attack  for  lift,  drag,  and  moment: 

<*re  (deg) 

maximum  peak  dynamic  stall  vortex  induced  lift  coefficient: 
AcSds 

maximum  peak  dynamic  stall  vortex  induced  drag  coefficient: 
A°dds 

maximum  peak  dynamic  stall  vortex  induced  moment  coefficient: 
ACmds 
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NLRTR 


KHLM.DA 

KFLMDA 

FXLHDA 

FYLMDA 

FMLMDA 

FACTWU 

KINTH 

KINTF 

KINTHB 

KINTHT 

KINTVT 

HINGE 

NGOLB 

NGCLT 

NONROT 

EPMODE 


I  exC"tO  r  k ^  for  hover  induced  velocity  (typically  l.l) 

factor  kf  for  forward  flight  induced  velocity  (typically  1.2) 

factor  f  for  linear  inflow  variation  in  forward  flight 
(typically  1.5) 

factor  f  for  linear  inflow  variation  in  forward  flight 
(typically  1.) 

factor  f  on  linear  inflow  variation  due  to  hub  moment 
(typically  1.) 

factor  introducing  lag  in  C„,  ,  and  used  to  calculate 

induced  velocity  (typically  .5)  X  ^ 

factor  for  hover  interference  velocity  at  other  rotor 
(k-21  or  k12) 

factor  for  forward  flight  interference  velocity  at  other 
rotor  (K21  or  K12) 

(linear  variation  between  KINTH  at  lv  =  O.O5  and  KINTF  at 
/•  =  0.10  is  used)  / 

factor  for  rotor- induced  interference  velocity  at 
wing-body,  Ky 

factor  for  rotor-induced  interference  velocity  at 
horizontal  tail, 

factor  for  rotor-induced  interference  velocity  at 
vertical  tail,  Ky 

(Ky,  Kjj,  Ky  equal  fraction  of  fully-developed  wake  times 
maximum  fraction  surface  in  wake) 

integer  parameter  specifying  blade  mode  type 
0  hinged 

1  cantilever 

2  articulated  (flap  and  lag  modes  only) 

number  of  collocation  functions  for  bending  mode  calculations 
(total  flap  and  lag,  alternating);  maximum  20 

number  of  collocation  functions  for  torsion  mode  calculations; 
maximum  10 

integer  parameters  NE  0  to  calculate  nonrotating  bending 
frequencies 

criterion  on  change  of  collective  pitch  to  update  blade 
modes,  (deg) 
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NLRTR 


MASST 

XIT 

MB  LA  EG 

EFLAP 

ELAG 

KFLAP 

KLAG 

RCPLS 

TSPRNG 

RCFL 

NOPB 

WTIN 

FTO 

FTC 

FTR 

KTO 

KTG 

KTR 

KPIN 


PHIPH 

PHIPL 

RPB 

RPH 

XPH 


tip  mass  (slug  or  kg);  the  tip  mass  can  also  be  included 
directly  in  the  section  mass  distribution 

offset  of  tip  mass  center  of  gravity  aft  of  elastic 
axis,  Xp/R 

blade  mass  (slug  or  kg);  if  LE  0.,  integral  of  section  mass 
used  (with  mass  included  at  r  =  0.  +o  account  for  the  hub  mass) 

flap  hinge  offset  ep/R  (extent  of  rigid  hub  for  cantilver  blade) 

lag  hinge  offset  ep/R  (extent  of  rigid  hub  for  cantilver  blade) 

flap  hinge  spring  (ft- lb/rad  or  m-N/rad) 

lag  hinge  spring  (ft-lb/rad  or  m-N/rad) 

hinge  spring  parameter,  ft. 

s 

hinge  spring  parameter ,  © 

SO 

(hinge  spring  pitch  angle  is  ©  =  ©  +  JR  ©„_) 

S  SO  S  ( j 

structural  coupling  parameter  JR  (effective  pitch  angle  JR® 
used  to  calculate  blade  bending  modes;  normally  JR=  1.) 

integer  parameter  specifying  twist  inboard  of  rpA:  EQ  1  for 
no  pitch  bearing 

integer  parameter  defining  control  system  stiffness  input: 

1  for  ,  2  for 

control  system  frequency  (per  rev,  at  normal  tip  speed  VTIPN) 
collective 

cyclic 

reactionless 

control  system  stiffness  (ft-lb/rad  or  m-N/rad) 
collective 

cyclic 

reactionless 

integer  parameter  defining  pitch/bending  coupling  input: 

1  for  input,  2  for  calculated  (negative  to  suppress  cosine 
factors  in  Kp^  and  Kp^,) 

root  geometry  to  calculate  pitch/bending  coupling  (KPIN  =  2  or  -2) 
pitch  horn  cant  angle,  <£>p^  (deg) 

pitch  link  cant  angle,  (deg) 

pitch  bearing  radial  location,  ipg/R 

pitch  horn  radial  location,  rp^/R 

pitch  horn  length,  xp^/R 


CO 


0 


atankp(nbm) 

DELJG 

RFA 

ZFA 

XFA 

CONE 

DROOP 

SWEEP 

FDROOP 

FSWEEP 


(de^  *** 

pnch/gl„hal  coupling  tan-%0  (deg) .  for  pitch  hor„ 

feathering  axis  radial  location,  r  /r 
gimbal  undersling,  Zp  /R  A 

torque  offset,  x^/r 
precone  angle  iFA{  (deg),  positive  up 

precone'®16  (deg)  at  S75  "  PMlttve  down  fro. 

sweep  angle  iFA}  (deg)  at  &?5  .  0,  positive  aft 
SomthS"L“1S  toCP  a"8le  S™4  <*>*).  positive  down 
feathering  axis  sweep  angle  SFAj  (deg),  positive  aft 


NLRTR 
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Namelist  NLWAKE 


FACTWN  factor  introducing  lag  in  bound  circulation  used  to  calculate 
induced  velocity 

OfVXVY  integer  parameters  EQ  0  to  suppress  x  and  y  components  of 
induced  velocity  calculated  at  the  rotors 


KNW 

KRW 

KFW 

KEW 


RRU 

FRU 

PRU 

FNW 

DVS 

DLS 

C0RE(5) 


0PC0RE(2) 


0PNWS(2) 


extent  of  near  wake,  Kj^ 

extent  of  rolling  up  wake, 

extent  of  far  wake  and  tip  vortices,  Kp^ 

extent  of  far  wake  and  tip  vortices  for  points  off  rotor 
disk, 

(age  4  =  KAH>  ;  all  K  GE  l) 

initial  radial  station  of  wake  rollup,  r^/R 

initial  tip  vortex  fraction  of  P  „  for  rollup,  fDtI 

max  nU 

extent  of  rollup  in  wake  age,  *RU  (deg) 

tip  vortex  fraction  of  ^  for  near  wake,  f^ 

sheet  edge  test  parameter  dvg}  LT  0.  to  suppress  test 

lifting  surface  correction  parameter  dls5  LT  0.  to  suppress 
correction 

vortex  core  radii  r^R 

( 1 )  tipvortices 

(2)  burst  tip  vortices 

(3;  tip  vortices  in  far  wake  off  rotor 

(4)  trailed  lines  (LT  0.  for  default  =  s/2) 

(5)  shed  lines  (LT  0.  for  default  =  t/2) 

integer  parameter  specifying  vortex  core  type:  0  for 
distributed  vorticity,  1  for  concentrated  vorticity 

(1)  tip  vortices 

(2)  inboard  wake 

integer  parameter  controlling  action  when  inflow  and 
circulation  points  coincide  in  near  wake  (  4  =  0)  an^ 
sheets  are  being  used:  0  to  use  two  sheets,  1  to  use  lines, 
2  to  use  single  sheet 

(1)  shod  wake 

(2)  trailed  wake 


LHW  number  of  spirals  of  far  wake  for  axisyrametric  case,  L^ 

OPHW  integer  parameter:  EQ  0  for  axisymmetric  wake  geometry 

OFRTS  Integer  parameter:  NE  0  to  include  rotation  matrices 

(Rrpg,  etc.)  in  influence  coefficients 
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NLWAKE 


WKM0DL(l3)  integer  parameter  defining  wake  model:  0  to  omit 

element,  1  for  line  segment  with  stepped  circulation 
distribution,  2  for  line  segment  with  linear  circulation 
distribution,  3  for  vortex  sheet  element 

(1)  tip  vortices  (stepped  line  or  linear  line) 

(2)  near  wake  shed  vorticity 

!3)  near  wake  trailed  vorticity 
h)  rolling  up  wake  shed  vorticity 
5)  rolling  up  wake  trailed  vorticity 

(6)  far  wake  shed  vorticity 

(7)  far  wake  trailed  vorticity 

(8)  far  wake  (off  rotor)  shed  vorticity 

(9)  far  wake  (off  rotor)  trailed  vorticity 

(io)  bound  vortices  (no  sheet  model) 

(11)  axisymmetrical  wake  axial  vorticity  (no  line  model) 

(12)  axisymmetrical  wake  shed  vorticity  (no  line  model) 

(13)  axisymmetrical  wake  ring  vorticity  (no  line  model) 


MRG 

NG(MRG) 


MRL 

NL(MRL) 

0PWKBP(3) 


number  of  circulation  points  for  near  wake;  LE  MRA 

circulation  points,  identified  by  aerodynamic  segment 

number:  ng^  for  i  =  1  to  MRG  (corresponding  r^  must  be 

between  r  ,/R  and  1.) 

root'  ' 

number  of  inflow  points;  LE  MRA 

points  at  which  the  induced  velocity  is  calculated, 
identified  by  aerodynamic  segment  number:  ng,  for 
i  =  1  to  MRL 

integer  parameter  controlling  blade  position  model  for  wake 
analysis 

(1)  EQ  0  to  suppress  inplane  motion 

(2)  EQ  0  to  suppress  all  harmonics  except  mean 

(3)  EQ  0  for  linear  from  r  =  r  /R  to  r  =  1. 


NLWAKE 


VELB 

DPHIB 

DBV 

QDEBUG 


KRWG 

OPRWG 


FWGT( 2) 

FWGSI(2) 

FWGSO(2) 

KWGT(4) 

KWGSl(4) 

KWGSO(4) 


core  burst  propagation  rate,  13 

core  burst  age  increment,  (deg) 

core  burst  test  parameter  d^v j  LT  0.  to  suppress  bursting 

velocity  criterion  for  debug  print:  print  if 
Iv’-k/r  \  >  Q.BEBUG 


Prescribed  wake  geometry 

extent  of  prescribed  wake  geometry,  K^y.,  (age  4*"  KA-V); 
maximum  144 

integer  parameter  defining  prescribed  wake  geometry  model 
1  from  K^f^X  *  K2=f2^  ’  input  K3»  input  K4 
option  #1,  without  interference  velocity  in  X 


2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


from  input  ,  Kg.  K^, 

Land grebe  prescribed  wake  geometry 

from  C,p 

from  P_ 

max 

from  X 

from  V.  without  interference 
Kocurek  and  Tangier  prescribed  wake  geometry 
from  Cm 


from  r 


max 


from  \ 

from  X  without  interference 


Factors  f.  and  f?  for  prescribed  wake  geometry 
xtip  vortex 

inside  shee+  edge 

outside  sheet  edge 

Constants  K.,  K_,  K~,  K^  for  prescribed  wake  geometry 
tip  vortex 

inside  sheet  edge 

outside  sheet  edge 


NLWAKE 


Free  wake  geometry 

KFWG  extent  of  free  wake  geometry  distortion  calculation,  _ 

(age  <£>  -  KA^);  suggest  (.4/k)MPSI;  maximum  96, 
multiple  MPSI 

OPFWG  integer  parameter  defining  free  wake  geometry  model 

1  Scully  free  wake  geometry 

2  option  #1,  without  interference  velocity 

ITERWG  number  of  wake  geometry  iterations;  suggest  2  or  3 


FACTWG 

RTWG(2) 

WGM0DL(2) 

CCREWG(4) 

MRVBWG 

LDMWG 

NDMWG(HPSI) 
DQWG( 2) 


factor  introducing  lag  in  distortion  calculation  to 
improve  convergence;  suggest  0,5 

radial  station  r/R  of  trailed  vorticity 

(1)  inside  sheet  edge 

(2)  outside  sheet  edge,  or  trailed  line  (suggest  .4) 


integer  parameter  defining  wake  model i  0  to  omit,  1  for 
line  segment,  2  for  sheet  element 

(1)  inboard  trailed  wake  elements 

(2)  shed  wake  elements 

vortex  core  radii  r  J'R 

(1)  tip  vortices 

(2)  burst  tip  vortices  (LE  0.  for  default  = 
unburst  value) 

(3)  inboard  trailed  lines  (LE  0.  for  default  = 
i(RTWG(2)  -  RTWG(l))  ) 

(4)  shed  lines  (LE  0.  for  default  =  0.4aV) 

number  of  wake  revolutions  used  below  point  where  induced 
velocity  is  being  calculated;  suggest  2 

integer  parameter  S'™,*  general  update  every  increment 

in  boundary  age;  suggest  lSO0/*^  " ! 


integer  parameter  nT_f(  'V .) :  boundary  update  every  n~. 
increment  in  age,  function  of  =  jAH>,  j  =  1  to  nFSI; 
suggest  90°/a>P  fore  and  aft,  and  450/A.vP  on  sides 


incremental  velocity  criteria;  suggest  0.04^  to  o.oe\i 

(1)  near  wake  elements  defined  by 
1  Aq  \  >  DQWG(l) 

(2)  integrate  bound  vortex  line  in  time  over 
if  |Aq  I  >  DQWG(2) 


NLWAKE 


IPWGDB(2)  integer  parameters  controlling  debug  level  3  print 
of  wake  geometry  distortion 

(1)  IPRi  print  distortion  before  general 
update  every  IPR  *  A*P ;  03  0  to  suppress 

(2)  INPSi  print  distortion  after  each 
iteration  every  IKPS  *  EQ  0  to 
suppress?  last  iteration  printed  in  full 

QWGDB  parameter  controlling  debug  level  3  print»  Induced  velocity 

contribution  of  wako  element  printed  if  |A‘q)  >•  QWGEB; 
suggest  0.5^  to  1.0^ 
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TITLE( 20) 
WEIGHT 

IXX 

IYY 

IZZ 

IXY 

IXZ 

IYZ 


title  for  airframe  and  drive  train  date  (30  characters) 
aircraft  gross  weight  including  rotors  (lb  or  kg) 

p  p 

aircraft  moments  of  inertia  including  rotors  (slug-ft  or  kg-m  ) 


I 

zz 

I 

xy 

"■yz 


TRATIO 

CONFIG 


ASHAFT(2) 

ACANT(2) 

ATILT 

HMAST 

DPS  1 21 

CANTHT 

CANTVT 


ratio  of  rotor  ti2  rotational  speed  to  rotor  #1  rotational 
speed,  (transmission  gear  ratio  rj^/rj^) 

integer  parameter  specifying  helicopter  configuration 
0  for  one  rotor 

1  for  single  main  rotor  and  tail  rotor  (rotor  #2 
is  the  tail  rotor) 

2  for  tandem  main  rotors  (rotor  it 2  is  the  rear  rotor) 

3  for  tilting  proprotor  aircraft  (rotor  #2  is  the 
left  rotor) 

shaft  angle  of  attack  (deg),  positive  rearward 

(1)  rotor  tt\ 

(2)  rotor  #2 

shaft  cant  angle  <t>p  (deg)  j  positive  to  right  for  main 
rotor;  positive  upward  for  tail  rotor;  positive  inward  in 
helicopter  mode  for  tilt  rotor 

(1)  rotor  #1 

(2)  rotor  #2 

nacelle  tilt  angle  ©<p  (deg),  for  tilting  proprotor 
configuration  only;  0.  for  airplane  mode,  90.  for 
helicopter  mode 

rotor  mast  length  from  pivot  to  hub  (ft  or  m),  for 
tilting  proprotor  configuration  only 

(deg);  rotor  #2  azimuth  angle  when  rotor  #1 
azimuth  angle  =  0;  must  be  0.  if  £L.y£2.^  j/  1. 

horizontal  tail  cant  angle  (deg),  positive  to  left 

vertical  tail  cant  angle  cj>VT  (deg),  positive  to  right 
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rscc 

BLCG 

WLCG 

FSR1 

3LR1 

WLR1 

FSR2 

BLR  2 

WLR2 

FStfE 

BLWB 

WLWB 

FSHT 

BLHT 

WLHT 

FSVT 

BLVT 

WLVT 

FSOFF 

BLOFF 

WLOFF 


location  (fuselage  station,  butt  line,  and  waterline)  of 
aircraft  components  relative  to  a  body  fixed  reference 
system  having  an  arbitrary  orientation  and  origin;  fuselage 
station  (FS)  posi  ive  aft,  butt  line  (BL)  positive  to  right, 
and  waterline  (WL)  positive  up  (ft  or  m) 

aircraft  center  of  gravity  location 


rotor  #1  hub  location  (right  nacelle  pivot  location  for 
tilting  proprotor  configuration) 


rotor  #2  hub  location 


wing-body  center  of  action 


horizontal  tail  center  of  action 


vertical  tail  center  of  action 


point  off  rotor  disk  (for  induced  velocity  calculation) 


CNTRLZ(li) 


control  inputs  (deg)  for  all  sticks  centered  (v  =  o)s 

e.-  a  6> 


ic  is 


ic  is 


*f  K 


rotor  #1  rotor  #2  aircraft 


o;>*> 


NLBCDY 

description  of  control  system  (for  TCpp);  K  parameters 
are  gains  (deg  per  stick  deflection),  parameters  are 

swashplate  azimuth  lead  angles  (deg) 

one  rotor,  single  main  rotor  and  tail  rotor,  tilting  proprotor 
configurations 

KOCFE  Kq,  collective  stick  to  collective  pitch 

KCGFE  K  ,  lateral  cyclic  stick  to  cyclic  or  differential 

collective  pitch 

KSCFE  K  ,  longitudinal  cyclic  stick  to  cyclic  pitch 

s 

KPCFE  K  ,  pedal  to  tail  rotor  collective  or  differential 

cfclic  pitch 

PGGFE  ,  lateral  cyclic  stick  to  cyclic  pitch  (one  rotor, 

or  single  main  rotor  and  tail  rotor  configurations) 

PScFE  •A.'V  ,  longitudinal  cyclic  stick  to  cyclic  pitch 

s 

PPGFE  ,  pedal  to  differential  cyclic  pitch  (tilting 

proprotor  configuration  only) 

tandem  main  rotor  configuration 

KFOCFE  KpQ,  collective  stick  to  front  collective  pitch 

KROGFE  Krq,  collective  stick  to  rear  collective  pitch 

KFUGFE  KpC,  xateral  cyclic  stick  to  front  cyclic  pitch 

KRCCFE  K^,  lateral  cyclic  stick  to  rear  cyclic  pitch 

KFSGFE  Kpg,  longitudinal  cyclic  stick  to  front  collective  pitch 

KRSCFE  Krs,  longitudinal  cyclic  stick  to  rear  collective  pitch 

KFPCFE  KpP,  pedal  to  front  cyclic  pitch 

KRPGFE  KRp,  pedal  to  rear  cyclic  pitch 

PFCCFE  lateral  cyclic  stick  to  front  cyclic  pitch 

PRGCFE  lateral  cyclic  stick  to  rear  cyclic  pitch 

PFPCFE  AH^pp,  pedal  to  front  cyclic  pitch 

PRPGFE  Af^p,  pedal  to  rear  cyclic  pitch 

aircraft  controls  (ail  configurations) 

KFCFE  Kp,  collective  stick  to  flaperon 

KTCFE  K,,  collective  stick  to  throttle 

t 

KACFE  Ka,  lateral  cyclic  stick  to  ailerons 

KEGFE  Ke,  longitudinal  cyclic  stick  to  elevator 

KRCFE  K  ,  pedal  to  rudder 
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NEM  number  of  airframe  modes  for  which  data  supplied; 

maximum  10 

QMASS(NEM)  generalized  mass  including  rotors  (slug  or  kg) 

QFREQ(NEM)  generalized  frequency  U)k  (Hz) 

QDAMP(NEM)  structural  damping  g& 

QDAMPA(NEM)  aerodynamic  damping  Fq a  =  b  (Qj^/f^V2)/  Mqsk/v) 

(ft  or  m  )  k  k  > 

QCNTRL(4,NEM)  control  derivatives  F_  *  “o  ^(Qt/ik/2)/i&  for 
ie>  &a»  ^  (ft  /rad^or  mVrad 7 

OOFSYM(NEM)  integer  vector  designating  type  of  mode:  GT  0  for 

symmetric,  LT  0  for  antisymmetric;  only  required  for 
flutter  analysis  with  OPSYMM  NE  0 

-A 

ZETARl(3»NEM)  linear  mode  shape  at  rotor  #1  hub  (ft/ft  or  m/m) 

ZETAR2(3,NEM)  linear  mode  shape  ^  r  t.  rotor  #2  hub  (ft/ft  or  m/m) 

GAMARl(3,NEM)  angular  mode  shape  1?^  at  rotor  # 1  hub  (rad/ft  or  rad/m) 

GAMAR2(3,NEM)  angular  mode  shape  at  rotor  //2  hub  (  :ad/ft  or  rad/m) 

pitch/mast-bending  coupling  (rad/ft  or  rad/m) 

KPMCl(NEM)  KmCj^  “  -  for  rotor 

KPMSl(HEM)  KHsk  -  -  bO^/bq^  ftc  rotor  #1 

KPMC2(NEM)  K^Qk  =  -  bG^/^q^  for  rotor  #2 

KPMS2(NEM)  KMSk  -  -  qSk  for  rotor  ttZ 


NLBODY 


Aircraft  aerodynamic  characteristics 


LFTAW 

Wing- body 

L*/q 

(ft2/rad  or  m2/rad) 

LFTFW 

Hf/q 

(ft2/rad  or  m2/rad) 

LFTDW 

LVq 

(ft2/rad  or  m2/rad) 

AMAXW 

X 

max 

(deg) 

IWB 

iWB 

(deg) 

DRGOW 

fWB  =  D(/q 

(ft2  or  m2) 

DRGVW 

f 

vert- 

(ft2  or  m2) 

DRGIW 

■"eV  “  (  V^/q)/  &L/q)2 

)_1  (ft2  or  m2) 

DRGFW 

D0ts/i 

(ft2/rad  or  m2/rad) 

DRGDW 

V* 

(ft2/rad  or  m2/rad) 

MOMOW 

v« 

(ft3  or  m3) 

MOMAW 

( Y t3/rad  or  m3/rad) 

MOMFVJ 

(ft3 /rad  or  m3/rad) 

MOMDW 

«JF/q 

(ft3 /rad.  or  m3/rad) 

SIDES 

Y@/q 

(ft2/rad  or  m2/rad) 

SIDEP 

V4 

(ft3/rad  or  m3/rad) 

SIDFIR 

vy'q 

(ft3/rad  or  m3/rad) 

ROLLS 

,Vq 

(ft3/rad  or  m3/rad) 

ROLLP 

™x/q 

(ft*/rad  or  n/Vrad) 

ROLLR 

™x/q 

(ftV^ad  or  m Vrad ) 

ROLLA 

(ft3/rad  or  m3/rad) 

YAWS 

%/q 

(ft3/rad  or  m3/rad) 

YAWP 

V^p/q 

(ftVrad  or  mVrad) 

YA  WR 

VNz/q 

(ft  /rad  or  m4/rad) 

yawa 

V5 

(ft3/rad  or  m3/rad) 

LFTAH 

Horizontal  tail 

Le</q 

(ft2/rad  or  m2/rad) 

LFTEH 

Wq 

(ft2/rad  or  m2/rad) 

AMAXH 

^max 

(deg) 

IHT 

iHT 

(d“g) 
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Vertical  tail 

LFTAV  L^q  (ft2/rad  or  m2/rad) 

LFTRV  L^/q  (ft2/rad  or  m2/rad) 

AMAXV  °<max  (de g) 

IVT  iVT  (deg) 

Airframe  interference 

FETAIL  -  (  ^6/  ^(lA)  )_1  (ft2  or  m2) 

LHTAIL  horizontal  tail  length  for  6.  (ft  or  m) 

HVTAIL  vertical  tail  height  h^^,  for  C",  positive  up  (ft  or  m) 

OPTINT  integer  parameter  controlling  airframe/ tail  aerodynamic 

interference!  EQ  0  to  suppress  (e  =  0  and  q-  **  0) 
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ENGPOS 


IENG 

KMAST1 

KMAST2 

KICS 

KENG 

GSE 

GSI 

KEDAMP 

THRTLC 


KPGOVE 

KPGOV1 

KPGOV2 

KIGOVE 

KIGOV1 

KIG0V2 


Engine  and  drive  train  parameters 

integer  parameter  specifying  drive  train  configuration: 

0  one  rotor 

1  asymmetric,  engine  by  rotor  #1 

2  asymmetric,  engine  by  rotor  #2 

3  symmetric 

2 

engine  rotational  inertia  r^  I„,  ^for  both^engines  if 
symmetric  configuration  (slug-Tt  or  kg-in  ) 

drive  train  spring  constants  (ft-lb/rad  or  n-N/rad) 
rotor  #1  shaft,  Kj^  or 

rotor  #2  shaft,  KM_ 

^2  2 

interconnect  shaft,  rj^K^.  or  rl  ^1 
engine  shaft,  r^K-, 

Hi  H* 

engine  shaft  structural  damping  g  ( degree  of  freedom) 

S  6 

interconnect  shaft  .  structural  damping  gg  ( degree  of 
freedom) 

engine  damping  factor  K-;  typically  1.0  for  turboshaft 
engines,  or  10.  for  induction  electric  motors 

&Pg/<^©>t  (dimensional),  for  both  engines  if  symmetric 
configuration j  if  the  throttle  variable  is  only  used 
for  the  governor,  just  the  products 
Kp  i>PE/i©t  -  -  W/AVS 

Ki 

must  be  correct  (P  =  .flpQp  =  -^gQg) 

governor  proportional  feedback  gains  (sec) 
to  throttle,  Kp  =  -  &YS 

to  rotor  #1  collective,  Kp  = 

to  rotor  #2  collective,  Kp  =  &§/  dVs 

governor  integral  feedback  gains 
to  throttle,  K^  =  -  ^tyg 

to  rotor  #1  collective,  KT  =  &B/ 

to  rotor  #2  collective,  KT  =  <}9/dH> 

JL  S 
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T1GOVE 

T1GOV1 

T1GOV2 

T2GOVE 

T2GOV1 

T2G0V2 


governor  time  lag  X,  =  2 t/cj  (sec) 
throttle  1  v  n  v  ' 

rotor  #1 

rotor  #2 

governor  time  lag  X,  =  l/to2  (sec2) 
throttle  ~  n  v  BC  > 

rotor  #1 

rotor  #2 


Namelist  NLLOA D 


Airframe  vibration 


MV  IB 

number  of  stations  for  airframe  vibration 
calculation  and  print;  maximum  10;  LE  0  to 
suppress 

airframe  location  for  vibration  calculation  (ft  or  m 

fsveb(mvib) 

fuselage  station 

blvib(mvib) 

butt  line 

wlvib(mvib) 

waterline 

_  y 

ZETAV(3,NEM,MVIB) 

linear  mode  shape  V  at  airframe  vibration 
stations  (ft/ft  or  m/m) 

MALOAD 

integer  parameter  controlling  print  of  motion  and 
aerodynamics!  EQ  0  to  suppress;  LT  0  for  only  plots 

MHLOAD 

integer  parameter  controlling  print  of  hub  and 
control  loads:  EQ  0  to  suppress 

MRLOAD 

number  of  radial  stations  for  blade  section  load 
calculation  and  print;  maximum  20;  LE  0  to  suppress 

rload(mrload) 

blade  radial  stations  r/R  for  section  loads 

MHARML 

number  of  harmonics  in  loads  analysis;  maximum  30; 

LT  0  for  no  harmonic  analysis;  suggest  about  MPSl/3 

NPOLAR 

integer  parameter  n  for  polar  plots:  symbol  printed 
every  n-th  step 

NWKGMP(4) 

integer  parameter  controlling  wake  geometry  printer 
plot;  EQ  0  to  suppress 
(i)  top  view 
(2j  side  view 

(3)  back  view 

(4)  axial  convection 

MWKGMP 

number  of  azimuth  stations  at  which  wake  geometry 
plotted;  maximum  8;  LE  0  for  no  plots 

JWKGMP(MWKGMP) 

azimuth  stations  at  which  wake  geometry  plotted 

(S>  =  j 
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NPLOT(75) 


integer  parameter  controlling  printer-plots  of  motion 
and  aerodynamics:  0  for  no  plot,  1  for  time  history 
plot,  2  for  polar  plot,  3  for  both  (only  time  history 
available  for  1-4  and  68-75) 


ill 

(3 

W 

(5) 

(6 

(7 

(8 

(9 

(10) 

(n) 

(12) 

13 

(14) 

(15) 

(16) 

(17) 

(18) 
(19 
(20) 

21 

(22) 

(23) 

(24) 


(27 

(28 

(29 

(30; 

(31 

32 

(33) 

(34 

35 

(36) 

(37) 
(38 
(39 
(40) 

41 

(42) 

m 

(44) 


bending  notion 
torsion  motion 
maximum  circulation 
\  off  rotor 


M 

A 

cx 

cd 

cm 

fd radial 


up 

u<p 

UR 

U 

© 

lag 

flap 

^eff’  }ift 
drag 

moment 

Meff,  lift 

drag 

moment 


interference 


UG 

VG 

1% 

D/c 

M/c 

Dr/c 

Fx/c 

Fr/c 

Fz/c 

Vc 


-  Cm /▼" 
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not  used 
not  used 
not  used 
Cp/<r 

cPi/% 

V/r 

cP0Ar 
L  * 

D  * 

M  * 

Jr 

Fx 

Fr  * 

F„ 


=  T 


rr 

not  used 
not  used 
not  used 


‘int 


* 

* 


rotating  frame  root  loads 
nonrotating  frame  hub  loads 
rotating  frame  root  loads 
nonrotating  frame  hub  loads 
section  loads,  shaft  axes 
section  loads,  principal  axes 
section  loads,  shaft  axes 
section  loads,  principal  axes 


* 

* 


* 


* 


♦dimensional  quantities 


for  polar  plots,  last  digit  of  integer  part  of 
(value/ increment)  is  printed,  if  it  is  a  multiple 
of  NPOLAR;  the  plot  increment  is  defined  as  follows 
.01  plots  27-35 

.1  plots  6,  8-16,  24-26,  36-51 

1.  plots  5,  7,  17-23,  52-61 

10.  plots  62-67 
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CT,  rs 


KFATIG 

SENDUR(18) 

CMAT(lR) 

EXMAT(lP) 


parameter  K  in  fatigue  damage  calculation;  suggest 
3  or  4 

endurance  limit  ST  (dimensional  force  or  moment) 

material  constant  G 

material  exponent  M 

rotating  frame  root  loads 

(1)  inplane  shear  f 

( 2 )  axial  shear  f  x 

(3)  vertical  shear  f 

(4)  flap  moment  m 

(5)  lag  moment  mz 

(6)  control  moment  mc 
nonrotating  frame  hub  loads 

(?)  drag  force  H 

(8)  side  force  Y 

(9)  thrust  T 

( 10)  roll  moment  M 

v 11 )  pitch  moment  H 

(12)  torque  Q 

section  loads  (principal  axes) 

(13)  chord  shear  f 

(14)  axial  shear  f 

(15)  normal  shear  i 

(16)  flatwise  moment  mx 

(l?)  edgewise  moment  m^ 

(18)  torsion  moment  m^ 

the  S-N  curve  is  approximated  by  N  =  C/ (S/Sg-  l)M 

use  SE  LT  0.  or  C  LT  0.  to  suppress  damage  fraction 
calculation;  use  M  EQ,  0.  to  suppress  equivalent 
peak-to-peak  load  calculation  as  well 


NLLCAD 
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MNOISE 

RANGE( MNOISE) 
ELVATN( MNOISE) 

AZMUTK( MNOISE) 

MHARMN(3) 

MTIMEN(3) 

AXS(MRA) 

0PN0IS(4) 


NLLOAD 


Far  field  rotational  noise 

number  of  microphones j  maximum  10;  LE  0  for  no 
noise  analysis 

microphone  range  relative  hub  (ft  or  m) 

microphone  elevation  relative  hub  (deg),  positive 
above  rotor  disk 

microphone  azimuth  relative  hub  (deg),  defined  as 
for  rotor  azimuth 

number  of  harmonics 

(1)  in  noise  calculation;  maximum  500 

(2)  in  aerodynamic  load  harmonic  analysis 
(suggest  MPSl/3) 

(3)  in  print  of  noise  (LE  0  for  no  print) 

number  of  time  steps  (LE  0  to  suppress) 

(1)  in  period  of  noise  calculation;  maximum  500 

(2)  increment  in  noise  print 

(3)  increment  in  noise  plot 

r\ 

blade  cross  section  area  A„s/c  at  aerodynamic 
segments,  for  thickness  noise  calculation( typically 
0.685  times  thickness  ratio) 

integer  parameter  controlling  noise  calculation: 

0  to  suppress,  1  for  impulsive  chordwise  loading, 

2  for  distributed  chordwise  loading 
(l^  lift  noise 

(2)  drag  noise 

(3)  radial  force  noise 

(4)  thickness  noise 
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CPF  LOW 

CPSYMM 

CPF  DAN 
NBLDFL 
MPSIFC 

NINTPC 

MP3 ICC 

DALPHA 

DMACH 

GPUSLD 

DELTA 

OPRINT 

OPGRND 

KASGE 

OPSAS 

KCSAS 

KSSAS 

TCSAS 

TSSAS 


integer  parameter  specifying  analysis  type:  LT  0  for 
constant  coefficient  approximation!  EQ  0  for  axial  flow; 

GT  0  for  periodic  coefficients 

integer  parameter:  NE  0  for  symmetric  and  antisymmetric 
analyses  (tilting  proprotor  configuration  only) 

integer  parameter:  EQ  0  to  suppress  flight  dynamics  analysis 

integer  parameter:  EQ  1  for  Independent  rotor  blade  analysis 

number  of  azimuth  steps  in  period  for  nonaxial  flow,  periodic 
coefficient  analysis  (OPFLCW  GT  0);  =  360/(N  K)  for 

old  number  of  blades,  =  ?20/(N,  ,M)  for  even  number  of 

blades 

integer  parameter  specifying  numerical  integration  option 
for  periodic  coefficient  analysis  (CPFLOW  GT  0):  1  for 

modified  trapezoidal  method,  2  for  Runge-Kutta  method 

number  of  azimuth  stations  (per  revolution)  in  evaluation 
of  average  coefficients  for  constant  coefficient  approximation 
(CPFLOW  LT  0);  =  360% 

angle  oi  attack  increment  Ac*  (deg)  for  calculation  of 
c^,  cd,  and  cm  derivatives  in  aerodynamic  coefficients 

Mach  number  increment  AJ’l/M  for  calculation  of  c^,  c^,  and 
c^  derivatives  in  aerodynamic  coefficients 

integer  parameter  controlling  use  of  unsteady  lift  and  moment 
in  flutter  analysis:  0  to  suppress;  1  to  include;  2  for 
zero  in  stall  (150><  Wl  <  165°) 

control  and  motion  increment  for  aircraft  stability  derivative 
calculation  (dimensionless) 

integer  parameter:  EQ  0  to  suppress  rotor/body  aerodynamic 
interference  in  flutter  analysis 

integer  parameter  controlling  ground  effect  analysis:  EQ  0 
for  out  of  ground  effect,  NE  0  for  in  ground  effect 

factor  for  antisymmetric  ground  effect:  0.  to  suppress,  1.0 
for  unstable  roll  moment  due  to  ground  effect  (tilting  proprotor 
configuration  only) 

integer  parameter  controlling  use  of  SAS:  EQ  0  to  suppress 

lateral  SAS  gain  Kc  =  -  (deg/deg) 

longitudinal  SAS  gain  K  =  /  Ci©tr  (deg/deg) 

lateral  SAS  lead  time  *cc  (sec) 

longitudinal  SAS  lead  time  T  (sec) 

s 
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NLFLUT 

0PT0RS(2)  integer  parameter t  EQ  0  for  rigid  pitch  model  (infinite 
control  system  stiffness,  no  pn  degree  of  freedom) 

(1)  rotor  #1 

(2)  rotor  #2 

D0F(30)  integer  vector  defining  degrees  of  freedom  for  flutter 

analysis;  0  if  not  used,  1  if  used,  2  if  quasistatic  variable; 
order » 

(0  L\\  c'.>  CA  C i\  I’O  CV*  O'! 

rotor  #1 
rotor  #2 

bending  pitch/ torsion  gimbal  rotor  inflow 

(15)  (9)  teeter  speed 

airframe  ftp'Vc:  ^  cVn,**V«u  ^ 

rigid  body  flexible  engine  governor 

body  (10)  speed 


C0N(26)  integer  vector  defining  control  variables,  0  if  not  used; 
order: 

rotor  #1  ©o  a,c  &■<,  ... 
rotor  # 2 

pitch  (8) 

airframe  <5^  Sr 

pilot  &■£ 

GUC(3)  integer  vector  defining  gust  components,  0  if  not  used; 

order:  Uq,  vq,  wq 

for  a  two-bladed  rotor,  is  replaced  by 

'  UVJ  •  X 

there  are  N,, ,  rotor  pitch  control  variables;  except  for  a 

two-bladed  rotor,  which  has  the  4  variables  Q.  ,  0.  , 

u  lc  Is  1 


AN  TYK.'4) 


NLFLUT 


nsysan 


NSTEP 

NF'Pb^ 

freq(nfreq) 


NBPLOT 


NXPLT 

NVPLT 

NAMEX1  ( NXFLT) 

NAMEVP( NVPLT) 

NI/PLT 

NFJPLT 

nfiplt 


MSPLT 


rtbr  L 

issLrssrst  ss™*  tesks  in  u-“r 

(1)  eigenarulysis 

(2)  transfer  function  printer-plot 
}//  tine  history  printer-plot 

V4)  rms  gust  response 

Eigenanalysis 

O.  f?rfeieenvalues'  1  f0r  “c™*™ 

cx&cuvectors,  iu  or  11  for  seros  is  well 
static  response  calculated  if  NE  0 

SSuLted^Tf ' 53  f°r  “hiCh  response 

iculated,  LE  0  to  suppress;  maximum  100 


vector  cf  frequencies  (per  rev) 


Transfer  function  printer-plot 

“n?s"e":sth0di  “  11  fr°"  “«“•>  “ 

number  of  degrees  of  freedom  to  be  plotted;  maximum  80 
number  of  controls,  to  be  plotted;  maximum  29 

n^rignorSfMe  *°  rlotted  (‘"consistent 

names'ignored)1"01  to  be  ( ^Intent 

frequency  stens  per  decade 
exponent  (base  10)  o  ‘  beginning  frequency 
exponent  (base  10)  of  end  frequency 

(maximum  NF  =  (NFIPLT  -  NFOPLT)  * NDPLT  + l  =  l$l) 

magnitude  plot  scale:  if  1  nin+  r.oi=.+< 

o  _i.x  ,  ..  ..  11  P-1-0^  relative  maximum  valuo> 

~  ^  leactxlve  if  3,  plot  relative  10, 


NLFLUT 


Time  history  printer-plot 


NTPLOT 

control  input  types  1  for  step,  2  for  impulse,  3  for 
cosine  impulse,  4  for  sine  doublet,  5  for  square  impulse, 

6  for  square  doublet 

PERPLT 

period  T  for  ir. pulse  or  doublet  (sec) 

DTPLT 

time  step  (sec) 

TMXPLT 

maximum  time  (sec);  maximum  NXPLT*NVPLT*TMXPLT/DTFLT =  7200 

NXPLT 

number  of  degrees  of  freedom  to  be  plotted;  maximum  80 

NVPLT 

number  of  controls  to  be  plotted;  maximum  29 

NAMEXP( NXPLT) 

vector  of  variable  names  to  be  plotted  (inconsistent 
names  ignored) 

NAMEVP( NVPLT) 

vector  of  control  names  to  be  plotted  (inconsistent 
names  ignored) 

Rms  gust  response 

LGUST(MG) 

real  vector  of  gust  correlation  lengths:  GT  0,,  dimensional 
length  L  ("cG  =  L/2V);  EQ  0.,  set  L  =  400.;  LT  0., 
magnitude  is  aimensionless  correlation  time  t:_  (frequency 

u>  =  n/nc:G) 

mgust(mg) 

real  vector  of  gust  component  relative  magnitudes 

MG  =  number  of  gust  components;  maximum  3 

namexa(macc) 

vector  of  names  of  degrees  of  freedom  for  which 
acceleration  calculated;  last  3  equal  ACCB  for  body 
axis  acceleration  (all  3  or  none)  (inconsistent  names 
ignored) 

freqa(macc) 

vector  of  acceleration  break  frequencies  (Hz);  2/rev 
used  if  LT  0.;  in  same  order  as  NAMEXA 

MACC 

number  of  accelerometers;  LE  0  for  none;  maximum  83 

FSACC 

location  of  point  at  which  body  axis  acceleration 
calculated  (ft  or  m) 

fuselage  station 

BLACC 

butt  line 

WLAGC 

waterline 

A 

Z:  fACG(3,NEH) 

linear  mode  shape  ^  at  point  where  body  axis  acceleration 
calculated 

NAMEXR(3) 

names  of  (2.  ,  and  S  in  state  vector;  assumed 

that  f^s*  and  follow  immediately  (inconsistent 

names  ignored) 
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Variable  names  for  linear  system  analysis 


Degrees  of  freedom 

11  p.c  (3,4  ... 

• 

bending 

1B15 

1T1 

• 

<n 

©)4  Bw/jl 

pitch/torsion 

• 

• 

1T9 

1BGC 

$<Zc. 

gimbal/ teeter 

1BGS 

PSIS 

rotor  speed 

1LU 

Vu 

inflow 

1LX 

V* 

1LY 

Vl) 

ZB1 

• 

« 

bending 

2B15 

2T1 

« 

©?  ©!?e* ...  9ri 

pitch/ torsion 

• 

« 

2T9 

2BGG 

gimbal,/  teeter 

2BGS 

(3gS 

PSII 

^x 

rotor  speed 

2LU 

Xu 

inflow 

2LX 

Xx 

2LY 

Vl) 

PHIF 

4*c 

rigid  body 

THTF 

PSIF 

XF 

x? 

YF 

ZF 
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rotor  #1 


y 

rotor  #2 


V 

airframe 


QF1 


flexible  body 


airframe 


NLFLUT 


QF10 

PSIE 

% 

TGOV 

1G0V 

^  S  eyvr, 

2G0V 

Control 

variables 

ICO 

1C1C 

1C1S 

1C4 

• 

90  ©lt. 

» 

1C8 

2C0 

2C1C 

2C1S 

2C4 

0O  •»* 

engine  speed 
governor 


rotor  #1 


V 

rotor  #2 


V 


2G8 

DELF  6c 


airframe 


NLFLUT 


J* 


Gust  components 
UG  uG 

VG  vG 
WG  wQ 


For  the  rotor  names, 
replaced  as  follows, 
CONFIG  =  0 
1 
2 
3 
3 


the  leading  character  (i  or  2)  is 
depending  on  the  helicopter  configuration 
blank  (left  justified) 

M  or  T 
F  or  R 

R  or  L  (OPSYMM  = 

S  or  A  (OPSYMM  j 


l] 


For  a  two  bladed  rotor,  BGC  is  replaced  by  BT 


For  first  order  degrees  of  freedom,  the  only  state  is  the 
velocity,  hence  it  is  the  velocity  that  will  be  plotted 
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Namelist  NLSTAB 


NPRNTP  integer  parameter  controlling  performance  print  during 

stability  derivative  calculation:  IE  0  to  suppress 

NPRNTL  integer  parameter  controlling  loads  print  during  stability 

derivative  calculation:  LE  0  to  suppress 

ITERS  number  of  wake  influence  coefficient/motion  and  forces 

iterations 

OPLMDA  integer  parameter  controlling  induced  velocity  calculation: 

if  0,  update  influence  coefficients  and  inflow;  if  1,  suppress 
influence  coefficient  update;  if  2,  suppress  inflow  update 
(and  influence  coefficient  update) 

DELTA  control  and  motion  increment  for  stability  derivative 

calculation  (dimensionless) 

B0F(?)  integer  vector  defining  degrees  of  freedom,  0  if  not  used; 

order:  4^,  ©p,  H'p,  Xp,  yp,  Zp, 

C0N(l6)  integer  vector  defining  control  variables,  0  if  not  used; 
order : 

rotor 

rotor  §2  9<>  ©it- 

airframe  S>z,  Er  ©t 

pilot  &©  St  Sp  &t 


GUS(3)  integer  vector  defining  gust  components,  0  if  not  used; 

order:  uG,  v  ,  wG 

CPPRNT(4)  integer  parameters  controlling  stability  derivative  print, 
EQ  0  to  suppress: 

(1)  rotor  coefficient  form,  dimensionless 

(2)  rotor  coefficient  form,  dimensional 

(3)  stability  derivative  form,  dimensionless 
(u)  stability  derivative  form,  dimensional 

KCSAS  lateral  SAS  gain,  K  =  -  2>&c/d<|>p  (deg/deg) 

KSSAS  longitudinal  SAS  gain,  K  =  bS> /d©,.  (deg/deg) 

5  SI’ 

TCSAS  lateral  SAS  lead  time  (sec) 

v 

TSSAS  longitudinal  SAS  lead  time  *Cg  (sec) 
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NLSTAB 


EQTYPE( 12) 


ANTYFE(5) 


NSYSAN 

NSTEP 

NFREQ 


integer  parameter  specifying  equations  to  be 
analyzed,  EQ  0  to  suppress 

with  ,  with  SAS 

(1)  complete 

(2)  symmetric 

(3)  antisymmetric 

with  ,  without  SAS 

(4)  complete 

(5)  symmetric 

(6)  antisymmetric 

without  ,  with  SAS 

(?)  complete 

(R)  symmetric 

(9)  antisymmetric 

without  ,  without  SAS 

(10)  complete 

(11)  symmetric 

(12)  antisymmetric 


integer  parameter  specifying  tasks  in  linear  system 
analysis,  EQ  0  to  suppress 

(1)  eigenanalysis 

(2)  transfer  function  printer- plot 

(3)  time  history  printer- plot 

(4)  rms  gust  response 

(5)  numerical  integration  of  transient 


Eigenanalysis 

calculation  control:  0  for  eigenvalues,  1  for  eigenvalues 
and  eigenvectors;  10  or  11  for  zeros  as  well 

static  response  calculated  if  NE  0 

number  of  frequencies  for  which  frequency  response 
calculated;  LE  0  to  suppress;  maximum  100 

vector  of  frequencies  (per  rev) 


T?T>T?r\  l  N 

a  nr  iiujVc^,  ) 


«f,r  9 


NBPLOT 

namexp(nxflt) 

namevp(nvplt) 

NXPLT 

NVPLT 

NDPLT 

NF'OPLT 

NF1PLT 

MSFLT 


NTPLO'. 

PERPLT 

DTPLT 

TMXPLT 

NXFLT 

NVPLT 

NAMEXP( NXPLT) 
NAMEVP( NVPLT) 


LGUST(MG) 


mgust(mg) 


NLSTAB 


Transfer  fanction  printer-plot 

calculation  method:  if  1,  from  matrices;  if  2, 
from  poles  and  zeros 

vector  of  variable  names  to  be  plotted  (inconsistent 
names  ignored) 

vector  of  control  names  to  be  plotted  ^inconsistent 
names  ignored) 

number  of  degrees  of  freedom  to  be  plotted;  maximum  ? 

number  of  controls  to  be  plotted;  maximum  19 

frequency  steps  per  decade 

exponent  (base  10)  of  beginning  frequency 

exponent  (base  10)  of  end  frequency 

(maximum  NF  = (NF1FLT-  NFOPLT) *  NDPLT +  1  -  15l) 

magnitude  plot  scale:  if  1,  plot  relative  maximum 
value;  if  2,  plot  relative  10**K;  if  3*  plot  relative  10. 

Time  history  printer-plot 

control  input  type:  1  for  step,  2  for  impulse,  3  for 
cosine  impulse,  4  for  sine  doublet,  5  for  square  impulse, 

6  for  square  doublet 

period  T  for  impulse  or  doublet  (sec) 

time  step  (sec) 

maximum  time  (sec);  maximum  NXPLT* NVPLT*TMXPLT/DTPLT  =  ?200 

number  of  degrees  of  freedom  to  be  plotted;  maximum  7 

number  of  controls  to  be  plotted;  maximum  19 

vector  of  variable  names  to  be  plotted  (inconsistent 
names  ignored) 

vector  of  control  names  to  be  plotted  (inconsistent 
names  ignored) 


Rms  gust  response 

real  vector  of  gust  correlation  lengths:  GT  0., 
dimensional  length  L  (Tq  =  L/2V);  EQ  0.,  set  L  =  400.; 
LT  0.,  magnitude  is  dimensionless  correlation  time 
(frequency  U  =  11/ *c^) 

real  vector  of  gust  component  relative  magnitudes 
MG  =  number  of  gust  components,  maximum  3 


-73- 


namexa(macc) 

vector  of  names  of  degrees  of  freedom  for  which 
acceleration  calculated;  last  3  equal  ACCB  for  body 
axis  acceleration  (all  3  or  none)  (inconsistent  names 
ignored) 

freqa(macc) 

vector  of  acceleration  break  frequencies  (Hz);  2/rev 
used  if  LT  0.;  same  order  as  NAMEXA 

MACC 

number  of  accelerometers ;  LE  0  for  none ;  maximum  10 

FSAGC 

location  of  point  at  which  body  axis  acceleration 
calculated  (ft  or  m) 

fuselage  station 

BLACC 

butt  line 

WLACC 

waterline 

Numerical  integration  of  transient 

TSTEP 

time  step  in  numerical  integration  (sec) 

TMAX 

maximum  time  in  numerical  integration  (sec) 

NPRNTT 

integer  parameter  n»  transient  print  every  n-th 
integration  step;  LE  0  to  suppress 

OPPLOT 

integer  parameter  controlling  printer  plot  of  body 
motion;  EQ  0  to  suppress 

DOFPLT( 21 ) 

integer  vector  designating  variables  to  be  plotted, 

EQ  0  if  not  plotted;  order: 

4*  Vp  ^  t  ^  if.  \  4  ©pVf  i* y  ^  Op  n  y  ‘ip  % 

OPTRAN 

see  namelist  NLTRAN 

CTIME 

CMAG(5) 

GTIME 

GMAG(3) 

GDIST(  ) 

VELG 

PSIG 

OPGUST(3) 

“I 

7 

NLSTAB 


Variable  names  for  linear  system 


Degrees  of  freedom 

PHIF 

THTl 

PSIF 

XF 

YF 

> 

ZF 

FSIS 

% 

Control  variables 

ICO 

iCIC 

1C1S 

2C0 

2C1C 

2C1S 

DELF 

DELE 

DELA 

DELR 

CT 

DELO 

DELC 

DELS 

i>. 

DSLP 

DELT 

&X 

Gust 

components 

UG 

UG 

VG 

VG 

WG 

WG 

analysis 
rigid  body 

rotor  speed 

rotor 

rotor  ->'2 

aircraft 

pilot 
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nlstab 

Of  the  rotor  control  names,  the  leading  character  (l  or  2) 
replaced  as  follows ,  depending  on  the  helicopter  configuration 
lumig  -  0  blank  (left  justified) 

1  H  or  T 

2  F  or  R 

3  R  or  L 

1°*  f.^st  °rder  degrees  of  freedom  the  only  state  is  the 
velocity;  hence  it  is  the  velocity  that  will  be  plotted 
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Namelist  NLTRAN 


N PR NTT 

NPRNTP 

NPRNTL 

NR3TRT 

TSTEP 

Tl'iAX 

ITERT 

CPLHDA 


DOF ( ? ) 


CPSAS 

NCSAS 

KSSAS 

TCSaS 

TSSAS 

OPPLOT 

D0FPLT(2l) 


integer  parameter  n:  transient/performance/loads  print 
every  n-th  integration  step;  LE  0  to  suppress 

integer  parameter  controlling  performance  print:  LE  0  to 
suppress 

integer  parameter  controlling  loads  print:  LE  0  to  suppress 

integer  parameter  n:  restart  file  written  only  every  n-th 
integration  step;  LE  0  to  suppress 

time  step  in  numerical  integration  (sec) 

maximum  time  in  numerical  integration  (sec) 

number  of  wake  influence  coefficients/motion  and  forces 
iterations 

integer  parameter  controlling  induced  velocity  calculation: 
if  0,  update  influence  coefficients  and  inflow;  if  1,  suppress 
influence  coefficient  update;  if  2,  suppress  inflow  update 
(and  influence  coefficient  update) 

integer  vector  defining  degrees  of  freedom  in  numerical 
integration;  EQ  0  to  suppress  acceleration; 
order:  <f>F,  ©j,  ,  H>F,  xp,  yF,  ^ 


integer  parameter  controlling  use  of  SAS:  EQ  0  to  suppress 
lateral  SAS  gain,  Kc  =  -  (deg/deg) 

longitudinal  SAS  gain,  Kg  =  Jb&s/  d©p  (deg/deg) 
lateral  SAS  lead  time  ^  (sec) 
longitudinal  SAS  lead  time 


Ts  (sec) 


integer  parameter  controlling  printer  plot  of  body  motion: 

EQ  0  to  suppress 

integer  vector  designating  variables  to  be  plotted;  EQ  0  for 
not  plotted;  order: 

4>f  c>f  Vp.  X  c  Md  'bp  Vp  Xp  vy  -ip  'Vj  Bp  •Jv  1  $ 
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NI.TR  AN 


GPTRAN 

CTIME 

CMAG(5) 


GTIME 

GMAG(3) 


GDIST(2) 

VELG 

P3IG 

0PGUST(3) 


Transient  gust  and  control 

integer  parameter  specifying  transient  option;  1  for 
control;  2  for  uniform  gust;  3  for  convected  gust 

period  T  for  control  (sec) 

control  magnitude  vPq  =  (  SQ  £>c  ^  &p  &t)  (deg) 

defines  cosine  control  transient  with  period  T 
and  magnitude  vpQ 

period  T  for  uniform  gust  (sec) 
v  T 

gust  magnitude  gQ  =  (u^  v^  w^)  (ft/sec  or  m/sec) 

defines  cosine  uniform  gust  transient  with 
period  T  and  magnitude  g 


lengths  for  convected  gust  (ft  or  m) 

(1)  wavelength  L 

(2)  starting  position  Lq 

gust  convection  velocity  V  (ft/sec  or  m/sec) 


azimuth  angle  of  convected  gust  wave  front  ^  (deg) 


integer  parameters  defining  convected  gust  model 
(l)  '  EQ  0  to  not  use  V 

(2;  rotor  #1:  0  for  gust  at  hub,  1  for  over  disk 

(3)  rotor  it 2:  0  for  gust  at  hub,  1  for  over  disk 


defines  cosine  convected  gust  transient  with 
wavelength  L  and  magnitude  gQ ;  for  Lo  =  R  the 
wave  starts  at  edge  of  rotor  disk,  for  L  =  0. 
the  wave  starts  at  hub  —  assuming  the  aircraft 
center  of  gravity  is  directly  below  the  hub; 
convected  at  rate  V  relative  to  moving  aircraft 
if  V  is  not  used,  It  rate  V  relative  to  fixed 
frame  if  V  is  used  ® 

o 


NLTRAN 


"'ransient  gust  and  control  subroutines 

The  subroutine  CONTRL  calculates  the  transient  control  time 
history5  C(t).  The  subroutine  GUSTU  calculates  the  uniform 
gust  time  history,  G(t).  The  subroutine  GUSTC  calculates 
the  convected  gust  wave  shape,  C(x  ).  The  subroutines 
presently  calculate  a  cosine-impulie  gust: 

CONTRL  C(t)  =  |(1  -  cos2Trt/T) 

GUSTU  G(t)  =|(1  -  cos  2tt  t/T  ) 

GUSTC  G(x  )  =  |(1  -  cos  2tt(x  -L  )/L  ) 

6  S  0 

Other  transients  may  be  use'’  by  replacing  these  subroutines 
as  required . 
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Namelist  Inputs  for  Old  Job  (Restart) 


Namelist  NLTRIM 

ANTYPE(3) 

OERSAD(iO) 

debug( 25) 

NPRNTI 

Namelist  NLFLUT 

ANTYPE(4) 

NSYSAN 

NAKSXR(3) 

Namelist  NLSTAB 

0PPRNT('4) 

KCSAS 

KSSAS 

TCSAS 

TSSAS 

EQTYPE^ 12) 

ANTYFE(5) 

NSYSAN 

( 

OPGUST(3) 

Namelist  NLTRAN 

NFRNTT 

NPRNTP 

NFRNTL 

NRSTRT 

TMAX 


7 .  NOTES  ON  PRINTED  OUTPUT 

This  section  presents  notes  on  the  printed  output  of  the  program, 
particularly  regarding  the  units  of  the  variables  appearing  in  the  output. 

Print  of  Performance  (Program  PERF) 

Operating  conditions 

a)  motion:  1st  number  dimensionless,  2nd  number  dimensional 

1)  velocity  =  ft/sec  or  m/sec 

2)  dynamic  pressure,  q  =  lb/ft^  or  N/in 

3)  weight,  Cy/V  =  lb  or  N 

4)  bod^  motion  =  deg/sec,  ft/sec  or  m/sec 

5)  z  =  ft/sec  or  m/sec 

6)  ^  =  rpm 

b)  body  orientation  and  controls  in  deg 
Circulation  convergences 

a)  tolerance,  CG/S  in  C^/v*  form 

b)  G/E  =  ratio  error  to  tolerance  (  ^  1,  if  converged) 

Motion  convergences 

a)  tolerance,  BETA  (etc)  in  deg 

b)  BETA/E  (etc)  =  ratio  error  to  tolerance  (<?  1.  if  converged) 
Airframe  performances  section  *1.2.6 

a)  aerodynamic  loads:  dimensional 

b)  components 

l)  angles  in  deg 

?.)  loads,  q  dimensional 

3)  induced  velocity,  tota?.  velocity  dimensionless 
Gust  velocity:  dimensionless 
System  powers 

a)  dimensional  (HP);  number  in  parentneses  is  percent  total  power 

b)  climb  power  =  V cW 


"O, 


«0 


System  efficiency  parameters! 

a)  gross  weight,  W  =  lb  or  N 

b)  drag-rotor  =  =  (P^+PQ)/Vj  D/q-rotor  =  D ^ii?'/2; 

L/D-rotor  =  W/D 

c)  drag-total  -  Dtotal  -  PtotalAi  D/q-total  •  &total/i^' 
L/D- total  -  »/Dtotal 

d)  figure  of  merit  =  H  -  1  -  r„cm.ldeal/Ptotal 


<* 


$  <• 


<*4> 

■  o 

o 

.&• 

l? 

S$' 

'<> 

W 

I;”/ 


Print  of  Rotor  Loads  (Program  LCADRl) 

Print  aerodynamics  (function  r  and  H^) 

a)  dimensionless  quantitie.  generally,  angles  in  degrees 

b)  induced  velocity  in  nonrotating  shaft  axes  (\x,  — \y»  •'Xj 

c)  interference  induced  velocity  is  that  due  to  other  rotor 

d)  gust  components  in  velocity  axes 

Force/c  (dimensionless)! 

'  mean  '  0 

L/C  -  1”2 


=  -§lr(c/c  )cP 


mean' 

=  -|U2(c/c  )  c, 

2  '  '  mean'  d 


=  4u2(c2/c  )c  =  M/c 
2  '  '  mean'  m  '  mean 

2/ 


=  ?\)C(c/c  )C{}  =  P _ ,/c 

2  '  '  mean'  aradial  1 


D/C 
M/C 
DR/C 
FZ/C 
FX/C 
«A/C 

fr/c 

FRT/C  =  F c 


=  L/c 
=  D/c 


mean 


mean 


radial'  mean 

.  OT/S  -  VcBean  -  d(y^)/dr 

-  V°- 


mean 

=  M  /c 
a'  mean 

-  F  /c 
r  mean 


mean 


'  so, 


'  'bi'. 
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Forces  (dimensional) 


L  =  section  lift 

(lb/ft  or  N/m) 

D  =  section  drag 

(lb/ft  or  N/m) 

M  =  section  pitch  moment 

(ft- lb/ft  or  m- 

DR  =  section  radial  drag 

(lb/ft  or  N/m) 

FZ  =  Fz  =  dT/dr 

(lb/ft  or  N/m) 

FX  =  F 

x 

(lb/ft  or  N/m) 

MA  =  N 

3, 

(ft- lb/ft  or  m 

FR  =  F 

r 

(lb/ft  or  N/m) 

FRT  =  Fr 

(lb/ft  or  N/m) 

Blade  section  power:  section  5*2.1 

CP/S  =  d(Cp/<v  )/dr 

P  =  section  power  (HP/ft  or  HP/m) 


> 

$ 


Print  During  Stability  Derivative  Calculation  (Program  STA3M) 

a)  increment:  1st  number  dimensionless,  2nd  number  dimensional 

b)  motion  and  controls:  1st  number  dimensionless,  2nd  number 
dimensional 

1)  angular  velocity  =  deg/sec 

2)  linear  velocity,  gust  velocity  =  ft/sec  or  m/sec 

3)  =  rpm 

4)  Zp  =  ft/sec^  or  m/sec^ 

5)  controls  =  deg 


cLlissd.  f*c2?css !  mcinsn'ts  and  heiress  in  ^ a  1*02?™ 
(rotor  #1  parameters,  body  axes);  torque  in  —  3  CL/Va 
form  (rotor  ft  1  parameters) 


Print  of  Stability  Derivatives  (Frogram  STABD) 

Options : 

* 

a)  rotor  coefficient  form,  M  X  =  #20/v-a 

b)  stability  derivative  form,  X  (acceleration) 

c)  dimensionless  or  dimensional 
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Dimensions: 


a) 


b) 


force  or  moment 

forces 

moments 

torq  ue 

M*X  form 

*mkO?/s 

iNibn2 

NLH2 

X  form 

A 

°2 

n 

b  2 

-n 

(FF) 

(FM) 

(FQ) 

subscripts 


acceleration  (z) 

-  rfR 

(FA) 

angular  velocity 

=  sn 

linear  velocity 

-  O-R 

(FV) 

controls 

-  5V.3 

gust  velocity 

=  HR 

(FV) 

Print  During  Flight  Dynamics  Numerical  Integration  (Program  STABP) 

a)  controls  in  deg 

b)  gust  velocity:  1st  number  dimensionless,  2nd  number 
dimensional 

c)  aircraft  motion:  1st  number  dimensionless,  2nd  number 
dimensional 

1)  displacement  =  deg,  ft  or  m 

2)  velocity  =  deg/sec,  ft/sec  or  m/sec 

3)  acceleration  =  deg/sec  ,  g 

4)  inertial  axes  =  deg/sec,  g 

Print-  Transient  Solution  (Program  TRANT) 
a''  controls  in  deg 

b)  gust  velocity  dimensional 

c)  aircraft  motion:  1st  number  dimensionless,  2nd  number 
dimensional 

1)  displacement  =  deg,  ft  or  m 

2)  velocity  =  deg/sec,  ft/sec  or  m/sec 

3)  acceleration  =  deg/sec  ,  g 

4)  inertial  axes  =  deg/sec,  g 
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d)  generalized  forces:  moments  and  forces  in  ^2C/V-a  form 
(rotor  it  1  parameters,  body  axes);  torque  in  —  tfCL/v-a 
form  (rotor  #1  parameters)  ^ 


8 .  UNITS 


The  program  will  work  with  English  or  metric  (Si)  units  for 
input  and  output.  Some  of  the  input  parameters  and  most  of  the  internal 
program  parameters  are  dimensionless  (based  on  the  rotor  radius,  the 
rotor  rotational  speed,  and  the  air  density).  The  units  for  input  and 
output  parameters  are  based  on  the  consistent  mass-length- time  system 
(foot-slug-second  or  meter-kilogram-second),  with  the  following 
exceptions : 

a)  The  aircraft  gross  weight  is  input  in  pounds  or  kilograms. 

b)  The  aircraft  velocity  is  input  in  knots  for  both  systems 

of  units  (alternatively  the  dimensionless  speed  can  be  input). 

c)  Power  is  output  in  horsepower  for  both  systems  of  uni+s. 

The  "dimensional"  output  for  angles  is  in  degrees;  the  "dimensionless" 
form  for  angles  is  in  radians. 
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9.  AIRFOIL  TABLE  PREPARATION 


This  section  describes  a  program  that  constructs  airfoil  table 
files  in  the  form  required  by  the  rotor  analysis.  The  program  will  also 
print  or  printer-plot  the  airfoil  data  in  the  file  being  created  or 
in  an  existing  file.  The  airfoil  tables  are  constructed  using  either 
analytical  expressions  or  an  airfoil  table  deck  (in  C81  format).  The 
subprogram  functions  and  namelist  input  labels  are  summarized  below. 

Subprogram 

Name 

MAINTB  Airfoil  table  preparation  (main  program) 

AEROT  Interpolate  airfoil  tables 

AEROPP  Printer-plot  airfoil  aerodynamic  characteristics 

Namelist 

Label 

NLTABL  fable  and  print/plot  data 

NLCHAR  Airfoil  characteristics  data 

The  structure  of  a  job  to  run  the  airfoil  table  preparation  program  is 
defined  below.  The  basic  structure  consists  of  the  following  steps;. 

1)  Airfoil  file  definition 

2)  Main  program  call 

3)  Title  card 

4)  Namelist  NLTABL 

5)  For  each  radial  station  (OPREAD  f  0),  either 

a)  Namelist  NLCHAR  (OPREAD  =  l) 

b)  Airfoil  table  card  deck  (OPREAD  =  2) 

Sample  jobs  are  presented  below. 

Create  airfoil  table  using  analytical  expressions. 

DDEF  FT40F001, .AIRFOIL 
CALL  MA1NPR0G 
title  card 

& NLTABL  table  data , NFAF=40,0PREAD=1 ,&END 
&NLCHAR  airfoil  characteristics  data,&END 
%END 
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Create  airfoil  table  using  C8l  format  airfoil  card  deck 

DDEF  FT40F001,, AIRFOIL 
CALL  MAINPROG 
title  card 

&NLTABL  table  data  ,NFAF=40,OFREAD=2,&END 
airfoil  card  deck 


Print  and  plot  airfoil  table  data 

DDEF  FT40F001,, AIRFOIL 
CALL  MAINPROG 
blank  card 

&NLTABL  output  data , NFAF=40 ,0FREAD=G  ,&END 

‘•'SEND 


The  following  pages  described  the  input  variables  and  data  for 
the  airfoil  table  preparation  program. 

First  Card 

TITLE(20)  title  (80  characters);  blank  card  for  OPREAD  EQ  0 

Namelist  NLTABL 


angle  of  attack  boundaries 

NAB  number  of  boundaries,  N  ;  maximum  20 

a 

NA(NAB)  indicies  at  boundaries,  n^ 

A(NAB)  at  boundaries  (deg,  -180°  to  180°) 

Mach  number  boundaries 

NMB  number  of  boundaries,  N  ;  maximum  20 

m 

NM(NMB)  indiciec  at  boundaries,  n^ 

M(NMB)  M  at  boundaries  (0,  to  1.) 

radial  segments 

NRB  number  of  segments,  N  ;  maximum  10 

R(NRB+l)  boundaries  of  segments  (R(l)=0,,  R(NRB+l)=l.) 

maximum  NAB*NMB*NRB  =  5000 
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0PPRNT(3) 


NMPRNT 

MPRNT(NMPRNT) 

NAPRNT 

AFkNT(NAPRNT) 

NFAF 


integer  parameter  controlling  output;  EQ  0  to 
suppress;  default  value  is  1 

(l)  interpolate  and  print 
f2)  interpolate  and  plot 
(3)  list  tables 

number  of  Mach  number  values  for  print  and  plot; 
maximum  10 

Mach  number  values  for  print  and  plot 

number  of  angle  of  attack  values  for  print;  maximum  60 

angle  of  attack  values  (deg) 

unit  number  for  airfoil  table  file  (default  4o) 


OPREAD  integer  parameter:  EQ  0  to  read  airfoil  table  and 

print  data  only;  EQ  1  to  create  airfoil  table  using 
analytical  expressions,  write  airfoil  file,  and 
print  data  (default);  EQ  2  to  create  airfoil  table 
using  G81  format  airfoil  card  deck,  write  airfoil 
file,  and  print  data 


Namelist  NLCHAR  (for  each  radial  station;  if  OPREAD  -  l) 


GLA 

a  =  c.  at  M  =  0  (per  rad)  (default  5.7) 

MDIV 

lift  divergence  Mach  number  Mdiv  (default 

OLMAX 

c*max  M  =  0  (default  1.2) 

FSTALL 

factor  f  for  Co  (default  0.5) 

s  •‘'max  v  ‘ 

MSTALL 

Mach  number  M  for  co  (default  0.4) 

s  *max 

GSTALL 

factor  gg  for  stall  c^  (default  1.2) 

HSTALL 

factor  h  for  stall  ce  (default  0.4) 

S  -Y 

CLF 

c^  for  stall  c^  (default  1.12) 

CMAG 

-mac  v •  / 

CMS 

cms  (default  -0.07) 

DELO 

SQ  (default  0.0084) 

DELI 

(default  -0.0102) 

DEL  2 

(default  0.384) 

DCDDM 

^Cd/^M  (default  0,65) 

MCRIT 

critical  Mach  number  at  £<.=  0  (default  0 

ACRIT 

critical  Mach  number  zero  at  «=  ^crit  1 

ALFD 

drag  stall  angle  (deg)  (default  10.) 

CDF 

cdf  for  stall  c^  (default  2.05) 

(default  33*) 


Airfoil  Card  Deck  (for  each  radial  station;  if  OPREAD  =  2) 


I .  Header 

a)  Card  1,  format  (30A 1,612) 

title,  30  alphanumeric  characters 

NML,  number  of  Mach  number  entries  in  ca  table 

NAL,  number  of  angle  of  attack  entries  in  c$  table 
NMD,  number  of  Mach  number  entries  in  Co  table 
NAD,  number  of  angle  of  attack  entries  in  eg  table 

NMM,  number  of  Mach  number  entries  in  cm  table 

NAM,  number  of  angle  of  attack  entries  in  c^  table 

II.  Lift  Coefficient  Table 

b)  Card  2,  format  (?X,9F7.0);  2  or  more  cards  if  NML  ^  1C 

Mach  numbers  M.  to  MV_.T 
1  NML 

c)  Card  Ja.t  format  (F7.0,9F7.0) 

angle  of  attack, 

lift  coefficients  cv  at  M  -  ft,  to  M.mT  or  Mn 
X  1  NML  9 

Card  3b,  format  (7X*9F7.0);  1  or  more  cards  if  NML  ^  10 
lift  coefficients  c^  at  M  =  M^  to  M^^ 

d)  repeat  card  3  for  &<.-  to 

III.  Drag  Coefficient  Table 

e-g)  format  as  for  lift  coefficient  table 

IV.  Moment  Coefficient  Table 

h-j)  format  as  for  lift  coefficient  table 

V.  Parameter  Limits 

a)  M^  =  0;  data  extrapolated  for  M  >  M^;  Mach  numbers  in 
sequential  order 

b)  =  -180°,  &<  =  180°;  angles  of  attack  in  sequential  order 

c)  NM  >  2,  NA  ^  2  for  lift,  drag,  and  moment 

d)  (NM+l)(NA+l)<sr  501  for  lift,  1101  for  drag,  576  for  moment 
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For  OPREAD  -  1,  the  program  calculates  representative  airfoil 
characteristics  using  the  following  expressions  (refer  to  the 
accompanying  figures). 

A)  Below  stall 


a/Vl-M2 

=»„-•  a(l-n)/((l-«dlv)yi:iij’v  ) 


M  <  M , 


II,.  <  M  <  M ,  +.1 

div  div 

H  «^H..  +.1 
div 


CS  = 


C  =  Cr, 


°d  "  So  +  &i“  +  S2'k2  + 


Acd  =  max  (0,  (M-Mc)  ) 

M  =  max  (0.  M  . .  (l-  ,.)) 


B)  Stall  angle 


co  =  co  min 

AS  ■'max 


■  cAs/c>. 


^1.  fl-")  +  ^ 


G)  Stalled  lift  (!©<1  >  c<  ) 

s 


c  n  =  sign(o<) 


f(gs^s-  l«xQcgs  +  (<oc|-ocs)hscas 

L  gsws  "  *<s 

max  (*s  c*s  ’  c*f  sin 


CX  =  cXf  s*n  2  « 


if  M>  45c 


*>*  ■>  w«-«. 


D)  Stalled  moment  (l«X|  >  o<  ) 

s' 

cn  =  •[  sign(cx.)  (60-  +  (lex)  -^).75C|nf 


60-  <x 


sign(c<)  (90-  1  «1)  .?5cmf  +  (loci  -6o)< 


30 


mf 


I**  l  <  60° 

M>60° 


cmf  ~  “  4cd(<*  =90)  -  -Kcd(cx«o<d)+cdf) 
E)  Stalled  drag  (  !<*!  >  04,  ) 

A  ‘ 


cd  =  cd^  =  ^d^  +  cdf  sin  /  |o^~  ^d  90 

90  -  tv. 


) 


E)  Reverse  flow  (|ot\>90) 


use  effective  angle  of  attack  and  account  for  moment  axis  shift 
c<e  *  c*  —  ^  sign<a< 

Cm  =  cm  +  (2  Cos~e)c*  +  (2  sincxe)cd 
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